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Norms and exclusion theorems 
By 
F. L. BAUER and C. T. FIKE* 


The purpose of this paper is to state and prove some theorems which relate 
the characteristic roots of two matrices. 


1. Notation. Let ||x|| denote the norm of a finite dimensional vector x. Let 
|| A|| denote any norm of a square matrix A satisfying 


||4 =|] S|IAll [ll]. 


(Such a matrix norm is said to be “‘consistent with” the vector norm.) Let lub(A) 
denote the least upper bound of A, i.e., 


lub (A) = min {«: ||A x|| Sa||x| for all x}. 


It is easy to see that! 
lub(A) <||A||. 


For the concepts and notation used here, see HOUSEHOLDER [3] and OsTROw- 
SKI [5]. 

2. Let x be a characteristic vector of A and let A be the corresponding cha- 
racteristic root, i.e., Ax=Ax. Then, if B is an arbitrary matrix either A is a 
characteristic root of B or 

x= (AI — B)1(A — B)x. 
In the second case, we get 
||*|| Slub[(A7 — B)+(A — B))||x|| 
and therefore 


14S lub[(AI — B)1(A— B)J Ss [|(AI — B)1(A-— B)||. 
Therefore we have 


Theorem I.. The set {z: ||(z7 — B)+(A — B)||=1 or det (zJ — B) =0} contains 
all characteristic roots of A. 


‘The multiplicative property of the least upper bound can be used to derive 
the weaker condition 


1<lub[(A — B)] -lub(A — B) S||(A — B)>|] -|4 — BI]. 


From this follows 





* Oak Ridge National Laboratory. Operated by Union Carbide Corporation for 
U.S. Atomic Energy Commission. 

1 lub (A) is a norm, too. Moreover, lub (A B) S lub(A) lub(B). Thus lub(A) isa 
multiplicative norm, although ||A|| may not be. 

Numer. Math, Bd, 2 11 
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Theorem II. The set {z:||(zJ—B)||7S ||A—B|| or det(zJ—B)=0} 
contains all characteristic roots of A. 
Interesting special cases can be obtained from Theorems I and II by taking 


n" 
particular norms. For example, take ||x|| = max |x;,| (or \|a|| = > |x1) and 
the corresponding least upper bound i=1 


n n 
lub(A) = max > |a,,| (or lub (A) = max 5) |4,s|). 
i R= * int 


By taking B to be the diagonal of A, we get in this case the classical theorems 
of GERSHGORIN [1]. HOUSEHOLDER [2] has already used norms to prove these 
theorems. More general vector norms can be made to yield the more complicated 
results of LENG [4]. 


3. In foregoing sections, no restrictions were made on either A or B. If we 
make the assumption that B is normalizable, i.e., similar to a diagonal matrix A 
we can obtain the further results described below. Let 


B=PAP", 
Then 
AI—B=P(AI—A)P 
and 


{<lub[A (AI —A)“] lub (P) lub (P~}) ub(A — B) <||G(I— B) >|] x( P)||A — BI], 


where 
x (P) =lub(P)-lub(P*). 


(Note that x(P)=>1, since for any vector z=-0 


z= P Pz implies ||z|| S lub(P) - lub(P*) - ||z|].) 
We have now 


Theorem III. The set {z: ||(27 —A)+||*Sx(P) ||A — B|| or det (zI — B)=0} 
contains all characteristic roots of A. 

This theorem provides a means of studying the effect on the characteristic 
roots of small changes in the elements of a matrix. In this context, A —B 
would probably be regarded as an error matrix. For applications the result 
leaves something to be desired, since the same bound is given for every cha- 
racteristic root. 


4. We now start from Theorem II again and make use of the spectral de- 
composition of a matrix B having simple characteristic roots: 


B= > A, E;. 
k=1 


Note that B is normalizable, in this case. The F, are the principal idempotents 
(Frobenius covariants) of B. The resolvent of B has the representation 


(zJ — B)? mm (z—A,)7E,. 
k= 


1 
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Applying Theorem IT, we get the set 
n 
M = {2:45 ||A— Bl] D|e—Al [| Eall or det (27 —B) =o}, 
k=1 


which is contained in the union of the m sets 
= {z:|z—A,| Sn||A — B]| -|| Zi|l}- 
This gives 
Theorem IV. The set M@ ¢ U Mi contains ali characteristic roots of A. 
If |A;—A,|>n ||A —B]| (WE, \|+||£;||), the intersection M;> M; is empty. 


M; has no intersection with U M, for ||A — B|| sufficiently small, e.g., 
j= 
j-i 





+ A—B in [Ai — Aa 
(*) || |<< m wei WEF IEall 


Then M;, contains exactly one characteristic root, since the set M; degenerates 
to the point A; as ||A — B||—>0. Now if 77 and z€ M,, then 


Jz—4| 2|4,— 4] —n||4 — Bll - |] Zill- 
Applying the condition (*) gives 








E;\| 
Je ASL Al El 
| E; || > Ed + WE 
a Fe [aj — Aj| 
and from the definition of M and Theorem IV 
1-4 — By) I Wl < E14 
|4— Bil sll 4 — Bll. 


Condition (*) makes the left side of this inequality positive. The result can 
now be expressed as 


Theorem V. If ||A — B|| satisfies condition (*), each of the sets 


¥ sae |4 — BI — SAM + Il 
M; =f: jz —A,| S|] E,|| - ; aaa - ae where at a 
contains exactly one characteristic root of A. This result implies that as 
||A — B||—>0, the circle around the characteristic root A; of B containing a 
characteristic root of A goes to zero with ||£,|| -||A —B||. Similar results, 
replacing |z—A,| by |z—A,|°, hold if A; belongs to an elementary divisor of 
degree v. Theorem V does not suffer from the deficiencies attributed to Theorem 
III, and it is important to appreciate the significance of this. The circles of 
Theorem V do not all have the same radii. 





5. Theorem III can be simplified if the vector norm belongs to a special 
class which includes the most used ones — the Hélder norms and diagonal 
transforms of them. If 

lub (D) = max | d 


iil 


11* 
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where D is a diagonal matrix, we call lub(A) axis-oriented. Then for axis 
oriented norms, we have 
Theorem IIIa. The set {z: min |z—A,| < x(P) ||A — B||} or equivalently the 


union of the circles 
{z: |z — 4;| S%(P) ||A — Bll} 


contains all characteristic roots of A. 


6. We now discuss the quantity ||£; il which occurs in Theorem IV and 
Theorem V. Note that E;= x; Vs , where x; is a right and y?/ is a left characteristic 
vector. Trace (E;,)=1 means yx,=1. Now if the vector norm is strictly homo- 
geneous, i.e., if || x|| = || - xl, then 


lub (E,) = min {a: | yf" x| - || x;|| S «|| «|]} 
= [lll Iolo, 


where ||y/!||p>=min {a: |yj/x| <a ||x|| for all x} i is the dual norm to the given 
vector norm. Define a cosine of two vectors a, b” by 


bra 
cos (a, bt —_—_—_————— 
0) = Tal 


Then 
lub (E;) = : 


cos (x;, yH?) 


is the reciprocal cosine of the left and right characteristic vectors. Moreover, 
lub(E;) gives for axis-oriented, strictly homogeneous norms a lower bound for 
x (P), since 


lub (E,) = ||¥;l] -|19#"llo = ||P ell lle P7 Ilo S *(P) [Leill lleF"ll. 
i = » 
max wa max lub(E,;) S%(P). 


Thus, if x(P)=1, lub(E;)=1 or x; and y# are dual vectors. We call a matrix B 

normal if this is true for all pairs of characteristic vectors, and we can consider 

max lub(E£;) a measure of nonnormality of a normalizable matrix, and x(P) an 
+ 


indication of this measure. 
For euclidean norm, the inequality of KANToRovi¢ [6] gives even the sharper 
result 
IP ell lle PI] fax (P) + Onin (P) 
lle Nefll =~ 20 nax (P) Omin (P) 


where Omax(P)=lub(P), Omin(P) = [{lub(P~)]+ the largest and smallest singular 
value, resp. 





Therefore, for euclidean norm 


1 -~ 1 1 
ee max lub(£;) S ; («(P) +a) ° 
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Tschebyscheff-Approximationen in kleinen Intervallen I* 
Approximation durch Polynome 
Von 
H. MAEHLy** und CH. WITZGALL 


Es sei 
B (n) 
die Schar der Polynome 
P,(%) = Pot Pixt--> +f, x” 


mit beliebigen reellen Koeffizienten #;. ,=-0 wird nicht gefordert. Ferner sei 
im Intervall [—b, b] eine positive stetige ,,Gewichtsfunktion’‘ v gegeben 


v(x) >0 fir —bSx*Sb. 


Dann besagt ein fundamentaler Satz von TSCHEBYSCHEFF ([3], [1]), daB es zu 
jeder in [—}), db] stetigen Funktion / genau ein Polynom P€$%(m) gibt, fiir das 
der Ausdruck 


:= Max 
|z|S6 








mal 


Uv Vv 











<a 


den kleinstméglichen Wert annimmt. Dieses Polynom P nennen wir die ,, 7scheby- 
scheff-A pproximierende’‘ oder auch kurz die ,,beste Approximierende‘‘ von } im 
Intervall [— 6, 6] mit der Gewichtsfunktion v. 


Ziel dieser Arbeit ist es Naherungswerte fiir die ,,Fehlerschranke*‘ 


Ja|:=|/P=t 


Vv 








der Tschebyscheff-Approximierenden P anzugeben. Ausgangspunkt hierzu ist 
die Tatsache, daB im Grenzfall unendlich kleiner Intervalle die ,,Fehlerkurve' 
zu einer regularen analytischen Funktion / 
thet | 

v 





6:= 
ein Tschebyscheff-Polynom (vgl. [2]) 
n+1 * 
(14 Baa() 


wird. Um dies zu beweisen, wird ein Satz iiber die Stetigkeit der Tschebyscheff- 
Approximationen hergeleitet, der auch an sich von Interesse ist. 





* Diese Arbeit wurde erméglicht durch die Unterstiitzung des Bureau of Ships, 
US Navy, und dessen Applied Mathematics Laboratory, David Taylor Model Basin, 
unter Contract No 2406(00) mit der Universitat Princeton. 

** Jetzt Syracuse University, Syracuse 10, New York (U.S.A.) 
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Die Ubertragung dieser Uberlegungen auf Tschebyscheff-Approximationen 
durch gebrochen rationale Funktionen ist nicht ohne Einschrankung méglich. 
Es treten Ausnahmefialle auf, in denen Analoga zu den in dieser Arbeit formu- 
lierten Satzen nicht gelten. Die Behandlung der gebrochen rationalen Approxi- 
mationen erfordert daher mehr Raum und soll in zwei Fortsetzungsarbeiten 
({4], [5]) durchgefiihrt werden. . 


§ 1. Die Stetigkeit der Tschebyscheff-Approximation durch Polynome 


Wir fragen, wie sich eine Abanderung der stetigen Funktion / und der Gewichts- 
funktion v auf die Tschebyscheff-Approximierende auswirkt. Da wir Approxi- 
mationstheorie im Sinne von TSCHEBYSCHEFF betreiben, legen wir als MaB fiir 
die Abweichung zweier Funktionen voneinander die Norm 

[If — gl] := Max| f(x) — g(*)| 
zugrunde. 

(1.1) Satz. Es gibt eine Konstante K, die nur von f,v,n und b abhdngt derart, 
dap fiir jede weitere stetige Funktion g und Gewichtsfunktion w mit 


1 1 
Heil |] —tett |] s 


||P — U]| SK(llf— gl] +A] [le — 2), 








gilt 


wober PE B(n) und UECB(n) Tschebyscheff-Approximierende von f und g mit den 
Gewichtsfunktionen v und w im Intervall [—b, b] sind. | A| ist die durch ||(P —f)/v|| 
erklarte Fehlerschranke. 

Beweis. Der in der Einleitung zitierte Satz von TSCHEBYSCHEFF ([3], [Z]) 
besagt nicht nur die Existenz und Eindeutigkeit der Tschebyscheff-Approxi- 
mierenden, sondern charakterisiert dariiber hinaus das beste Polynom durch 
folgende geometrische Eigenschaft der Fehlerkurve 6-P—): &£s gibt min- 
destens n+ 2 aufeinanderfolgende Stellen . 


—bS% << CH SD, 
an denen die Fehlerkurve den Betrag 
na 
v 
mit alternierendem Vorzeichen annimmt. Falls diese ,,kritischen Stellen‘’ x; nicht 


eindeutig bestimmt sind, wollen wir die Auswahl stets so treffen, daB x, die 
kleinste Abszisse ist, fiir welche 








gilt. Wir definieren dann 
(1.2) A=A(t,v) := P(%o) — f(%o) 





und haben 


(1.3) (— 1)* (P(x) — f(x) =Av(x), t=0,...,.m+4. 
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Analog definieren wir yo und yu: 





U ae 
w=pm(g,w) := ) af . 
Wir wollen jedoch aus Bequemlichkeit fiir das folgende 
A=0 


voraussetzen, was ohne Einschrankung méglich ist, da man ohne weiteres — /, — g 
statt /, g betrachten kann. Uber das Vorzeichen von yu kénnen wir dann aller- 
dings nicht mehr verfiigen. Ist A=0, so kénnen die kritischen Stellen willkiirlich 
gewahlt werden. 


Weil U Tschebyscheff-Approximierende von g ist, gilt 











l= eS Fo 
Wegen 
a=f — FHL + 1g) + (4 (v—w)) 
hat man . 
(1.4) Ju] sa+||2| (tell + ally — a). 


An den kritischen Stellen x; von f gilt nun 


(— 4)" (P(x) — #(%)) = el 
(— 1)' (U(x) — g(%)) S| z| @ 
Subtraktion liefert 
(— 1)* (P(x) — U(x)) 2 — lf — gl| + 40(x) —|a| w 
Daraus folgt unter Verwendung von (1.4) 
(—4)* (Plas) — U(x) = — (1+ lool] |S I) (le — ell + Alle — w)- 


Wegen unserer Voraussetzung in (1.1) gilt jedoch 


Heol || S|] s+ dole 








also schlieBlich 
(1.5) (—4)* (P(x) — U(x) = — ( 


Wir bendtigen nun folgenden 





|) (lf —gll + Alle — wI)). 





(1.6) Hilfssatz iiber Polynome. Zu je n+ 2 Punkten 
—bS%y< Hy << % 4,50 
gibt es eine Konstante H so, daB fiir jedes Polynom DE B(n) mit der Eigenschaft 
(1.7) (—1)'D(x)=>—a, aZO0, 1t=0,...,"4+1 


(s. Fig. 1) gilt 
||| < Ha. 
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Beweis. Fiir a=0 ist nichts zu zeigen; denn ein Polynom D€$(m), das an 
n-+ 2 aufeinanderfolgenden Stellen x; den Ungleichungen 


(1.8) (—1)'D(x) 20, 1t=0,...,m+1 


geniigt, verschwindet trivialerweise identisch (D kann nicht (m+ 1)-mal das 
Vorzeichen wechseln; deshalb muB D fiir mindestens ein x; verschwinden. Das 
Polynom D*:=(x—x,)4D erfillt jedoch 
wieder Ungleichungen vom Typ (1.8), usw.). 

Ist a>0, so betrachten wir die Polynome 


E*C®R(n), k=0,...,.n+4, ,* 
$B (n) n aan 


die durch die Gleichungen 





(1.9) (—1)'E*(x,) =—1, fir i+k 


eindeutig festgelegt sind. Die Polynome a E* 
, begrenzen‘‘ die Gesamtheit der Polynome D 


im folgenden Sinn *. Fig. 1. Polynom, (1.7) geniigend 


(1.10) Min aE*(x)<D(x)< Max aE*(x) fiir jedes x 


und alle Polynome D€ $(n) mit der Eigenschaft (1.7). 


Wir beschriénken uns darauf D<MaxaE* zu zeigen. Wir betrachten die 
Bereiche 
(~ oo, %4), (1, X3), (%3, Xs), a 


und schlieBlich als letzten 
(X,41, 0) falls m gerade 
(x,, 0) falls m ungerade. 
Diese Intervalle iiberdecken die reelle Achse mit Ausnahme der Punkte 
Bas Me, Bos 60 
An diesen Punkten ist jedoch DS Maxa F* trivial. 


Ist nun % € (%g;-1, %2;41), So betrachteu wir die Werte des Polynoms 


, G:=D—akE* 
an den Stellen 


x; fiir «+27, 
em ‘ae ; Vi < Vi41- 
x fir +=2), 
Es gilt fiir 727 


(— 1)*G(y,) = (—1)' G(x) =a + (— 1)*D(x,) S0. 





* Die Polynome D mit der Eigenschaft (1.7) bilden einen abgeschlossenen kon- 
vexen Kérper im linearen Raum $(m), der als Ecken gerade die Polynome @E* hat 
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Ist nun (—1)*/G(¥2;)=G(x)=0, so besitzt G die Eigenschaft (1.8) und ver- 
schwindet identisch; d.h. G(x) >0 ist unméglich. Also 
D(x) Sa@E*i(x) fiir %€ (%_;-1, Xej41)- 
Entsprechend zeigt man, daB 


D(x) Sa E°(x) fiir x€(— 9, x,) 
D(x) Sa E(x) fiir x€(*,4,, 0) falls m gerade 
D(x) SaE"**(x) fir x€(x,, 00) falls » ungerade. 
Aus (1.10) folgt 
||| < Max ||E]-a. 


Da H:=Max||E*|| nur von den Werten x; abhiangt, ist unser Hilfsatz (1.6) 
bewiesen. * 


Nun gilt nach (1.5) 
(— 1) D(a) := (— 1)! (P(x) — U(x) 2 — (2+ [lo 
also nach Hilfsatz (4.6) 








4) (l¢—ell+4Jo—wI)), 


|P— Us (2 +H] 








) (lel +All» — a), 


a 
v 





was den Beweis von Satz (1.1) vollendet. 
Eine direkte Folge von Satz (1.1) ist 
(1.11) Satz: Aus 


lim g=/, lim w,=v  gleichmaBig in [— }, b] 


folgt 
lim 6(g,, w,) = 6(f,v) gleichmapig in [— 6, 6). 


Aus Satz (1.11) laBt sich 
(1.12) Lim | A(g,, «,)| = lim |] (gs, w)|| = |], »)|| = [2 9) 


entnehmen. Ohne Betragstriche ist jedoch (1.12) unter Umstanden falsch. Das 
liegt an unserer Festlegung (1.2) des Vorzeichens von A mit Hilfe der kleinsten 
kritischen Stelle x). Es kann namlich vorkommen, daB die kleinsten kritischen 
Stellen der g, nicht gegen die kleinste kritische Stelle von / konvergieren, und 
daB somit das Vorzeichen von A ,,umklappt‘‘. Es bedarf einer Zusatzbedingung, 
um die Konvergenz von A(g;, w,) im Fall A(f, v) -=0 zu sichern. 

(1.13) Satz. Sind die kritischen Stellen von f beztiglich der Gewichtsfunktion v 
eindeutig bestimmt, und ist 


lim g=f, lim w,=v  gleichmafig in [— b, bj, 


so gilt 


lim A (8¢, 1) =A(f, v). 
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Beweis. Die Eindeutigkeit der kritischen Stellen erfordert 


(4.14) A(}, v) +0. 
Wir fiihren folgende Abkiirzungen ein: 


6,:= 0(g;,%),  69:= (7,2) 

Ayi=4(8:,™%), Ag: =Alf,v). 
Weiterhin sei 

—bsxn<ai<--<wi1sb 


ein Satz von kritischen Stellen von g, beziiglich w,, wobei, — wie in (41.2) ver- 
abredet —, xj die kleinstmégliche kritische Stelle ist. Ist {t};., 5.0 eine 
Nullfolge, so besitzen die Folgen {x{};_, .. konvergente Teilfolgen. Wir wahlen 
fiir jedes k eine solche konvergente Teilfolge 


{xji*},- a > Xp , 
wobei {s,,};~1,9,, eine Teilfolge der Nullfolge {t,};_ 1,5... ist. 
Die Grenzwerte x, erfiillen die Ungleichungen 
— bS%S sath S 3441 & 9. 


Aus der Stetigkeit von f und der Tatsache, daB die Grenzwerte / und v gleich- 
maBig erreicht werden, folgt leicht 


6o (Xx) = lim 6,,, (x4!*) = lim (— 1)* A,,, = (— 1)* lim A 


. YT See? 
p00 tk tk 


Damit ist die Existenz von lim A,,, sichergestellt, und wegen 
t—> co 


Stk 
| Jim 2,,, | = lim | Ass,| = | Ao] + 0. 


Das bedeutet aber gerade, daB 6, an den Stellen x, die Fehlerschranke mit alter- 
nierendem Vorzeichen erreicht. Die Punkte x, sind also ein Satz von kritischen 
Stellen von /. Nach Voraussetzung sind aber die kritischen Stellen von / ein- 
deutig bestimmt ; jede Folge {x‘{} kann also in Wahrheit nur einen Haufungspunkt 
besitzen, namlich die k-te kritische Stellle x, von /, muB also gegen x, konver- 
gieren. Es gilt nun speziell 

lim A= lim 6; (%) = 40(%0) = Ao, 
was zu beweisen war. 


§ 2. Grenzverhalten der Fehlerkurve bei Approximation durch Polynome 
Es seien nun / und die Gewichtsfunktion v in einer Umgebung des Nullpunkts 
vegulare analytische Funktionen 
f(%) = Sq + 81 ¥ + Sp %? +--+. 
v(x) =1+4,x4+4,x7+---. 
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Wir wollen zeigen, daB bei Tschebyscheff-Approximation durch Polynome im 
Grenzfall b->0, die ,,normierte Fehlerkurve‘‘ 


6(f,v) 
A(f,v) 





gegen das Tschebyscheff-Polynom 


(— 4)*** Tas (Z} 


strebt. Wir miissen diese Behauptung genauer formulieren. Mit der neuen Ver- 
anderlichen 
x 
71" = 
schreiben wir 
F,(y):= (by) = Sy +5, by + 520? y?+--- 
Vi(y) = v(by) =1+hby +h byt. 


Anstatt / und v im Intervall [— 0,0] betrachten wir F, und V, im Intervall 
[—1, +1] als Funktionen von y. Es sei 


6,:= 4(B,, V,) 
die Fehlerkurve zur Tschebyscheff-Approximierenden 
PE ®(n) 
von F, mit der Gewichtsfunktion V, und 
Ay: = A(K, Vs) 
die zugehérige Fehlerschranke. Wir behaupten dann 
(2.1) Satz: Ist s, ,;=——L-—— f+") (0) £0, so gilt 





(n+ 1)! 


lim +* =(—1)"7T,,  gleichmapig in [—1, +4]. 


Beweis. Wir setzen 
(2.2) F,(y) := so +s,by+---+5, "9" 
und betrachten die Funktion 


: 1 iG 
9° = —— saat (Fe — 5), 


wobei s,,, nach Voraussetzung von Null verschieden ist. g, entsteht somit 
aus #,, indem man +1 Glieder der Potenzreihe streicht und den ersten 
Koeffizienten der Restreihe auf eins normiert. Es gilt also 


I= Pr + Sab Yr + -.. 
und offensichtlich 
lim g»(y) = y"*?_ gleichmaBig in [— 1, + 4] 


lim V, (vy) =1 gleichmaBig in [— 1, + 41]. 
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Wir kénnen daher Satz (1.11) auf y”*?, 4 in der Rolle von /, v und gy, V, in der 
Rolle von g, w anwenden. Die Fehlerkurve, die sich bei Tschebyscheff-Approxi- 
mation von yt! durch ein Polynom n-ten Grades mit der Gewichtsfunktion 1 
ergibt, ist aber bekanntlich (vgl. [2]) das (n+ 1)-te Tschebyscheff-Polynom: 


d(y"*4, 1) = — $-Thsa(y). 
Daher liefert Satz (4.11) 
(2.3) lim 6 (gp, V,) =— 4 Tra (vy) gleichmaig in [— 1, + 4]. 


Wir haben nun den Zusammenhang zwischen der Fehlerkurve 6(g,, V,) und der 
Fehlerkurve 6,=6(F,, V,), die uns eigentlich interessiert, aufzuweisen. 


Es gilt auf Grund der Definition von 9, 











R-F, _ 
A. om f,—F, = F, — Fy— Sn410"t! m on pnt Sn410"t* he 
a 5 aa 0, |" 


Der Ausdruck in eckigen Klammern beschreibt die Fehlerkurve fiir die Approxi- 
mation von g, durch das Polynom gee 8 (x) mit der Gewichtsfunktion V,. 
Diese Fehlerkurve ist aber proportional zu 6,, hat also die geometrischen Eigen- 
schaften, die nach dem Satz von TSCHEBYSCHEFF die beste Approximation 


charakterisieren. Daher ist x et Tschebyscheff-Approximierende von g, 


und es gilt 
(2.4) bp = 8441" ** 3 (M5, Vi). 
Aus (2.3) und (2.4) folgt 


" T,,,, gleichmaBig in [— 1, +1]. 


. , = 
(2.5) lim ——Srr = jim 4(%, Vo) = 2 


Die kritischen Stellen von y"** sind eindeutig bestimmt, und zwar ist — 1 die 
kleinste kritische Stelle, so daB nach unserer Festsetzung (1.2) gilt 








" 1 wad 

A(y"*4, 4) = — 3 Traa(— 1) =P 
Satz (1.13) liefert nun 

hy mt+1 4) (—1)" 

lim a oo ia == lim A(p,,V,) = A(y"™", 1) = 
oder 
; Ay — (—1)"Sn41 

(2.6) lim prtti — 3" 


Damit sind wir aber fertig; denn (2.6) ergibt zusammen mit (2.5) 


lim % =(—1)"1 7... gleichmaBig in [— 1, + 4], 


b->0 Ay 


was zu beweisen war. 
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Als Nebenresultat haben wir (2.6) gewonnen. Es liefert gerade den in der 
Einleitung angekiindigten Naherungswert fiir 4,. (2.6) ist auch dann richtig, 
wenn S,,=0 ist. Fiir das in (2.2) erklarte Polynom F, gilt ja wegen der Opti- 
malitat von A, _ 

7 = 
lim ets im per| “yy 

Wir haben also 

(2.7) Satz. Sind f,v im Nullpunkt reguldre reelle analytische Funktionen 
und ist 








“b| = |s,41|=0. 


A, = A, (7, 2) 


die Fehlerschranke fiir die Tschebyscheff-Approximation von f durch ein Polynom 
P,€8(n) mit der Gewichtsfunktion v im Intervall [—b, +6], so gilt 





: —1)" Sp 1 oo 1 n 
few mr =t u . “i, ge f(0). 
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Uber einen Interpolations-Algorithmus 
und gewisse andere Formeln, die in der Theorie 
der Interpolation durch rationale Funktionen bestehen* 
Von 


PETER WYNN 


1. Einleitung 
1.1 Die Legendreschen Interpolations-Polynome 
Das klassische Problem, dasjenige Polynom n-ten Grades in der Variablen x 


zu konstruieren, das die Werte /f),/,,...,/, annimmt, wenn x die Werte %p, 
X,+-.,%, annimmt, wird auf folgende Weise gelést [1], [2]. Die Folgen divi- 
dierter Differenzen /,,(%9, %1,---, Xm—1» %;) M=0,1..., s=m, m+1,..., ent- 


stehen durch sukzessive Anwendung der Beziehungen 


fo(x.) =f, (s=0,1,...) 


fons (Bee ps +++ Syyn Hy) = Leas Sie =~ Sm—ae Fe) = fel os Sy ---» Sma tm) 4.4.4) 
X%3— Xm 





(m =0,1,...; s=m+1,m-+ 2,...). 
Das gewiinschte Polynom ist sodann 


fe,1,..,n-1,0(%|%e, xy; si atid Xn-1> Xn) 


=fot(*—%o) fi (%o, %1) +++ + (%— Xo) (%— 4%) © (¥—H% yp) fa (HO X15 +++ Xp) 


1.2 Die Aitkensche Methode 
Ein einfaches Schema zur Konstruktion der Polynomfolgen 
be,2....,9,0(%| See Bar +>» Spo B) (7 =0,1,...; $7 + 4,7 + 2,...) 


ist von AITKEN [3] gegeben worden, es wird jetzt in weitem Umfang zur prak- 
tischen Interpolation angewandt. Die Polynome werden in dem Schema 


%—%y fo (%| Xo) 
%—% A(%|%) fo (%|%0, m4) 
%—%_ fy(%|%2) fo,2(%|%0,%2) fo,r,2(%|%o, %1, %2) 


i fs(%| 3) fo,3(*| *0, 3) fo,1,3(%| Xo, %1, %3) fo,1,2,3(%| %o, %1 Xe, Xs) 





* Dissertation zur Erlangung des Doktorgrades bei der Naturwissenschaftlichen 
Fakultat der Johannes Gutenberg-Universitat zu Mainz. D 77. Der Verfasser méchte 
gerne Herrn Prof. Dr. F. L. Baukr fiir seine Ermutigung in dieser und anderen Unter- 
nehmungen herzlich danken. 
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angeordnet und sind durch die Beziehungen 


ae) | ee were 
__ (¥#— 45) fot, ...,7—1,5(¥ | 40s Mao ++» %e—1» Xs) — (¥—%s) fo,1,..., 1-1, 7 (¥ | Xo» Hay «++» Hp—1» Fr) (1.2.4) 
4%, — Xs 


(7 =0,1,...5 S=r+1,7+2,...) 





konstruierbar. 
1.3 Die Nevillesche Methode 
Die ,,zentralen’ Polynome fy mit,...,m+s(%|%m> Xm+i +++» Xm4s) Sind in ahn- 
licher Weise durch eine Methode, die von NEVILLE [4] gegeben worden ist, 
konstruierbar. Die Polynome werden in dem Schema 


fa (%) 

fo,1(*) 
h(®) hi,2(%) ae fo,1,2,3(*) 
fa(x) °° , 


angeordnet und entstehen durch die Beziehungen 


benn4',.... +043) 





> (¥ — Xm) fata, m+, .., mtr+1(¥%) — (¥ — ¥m4e+1) hm, m+1, ..., m+ (*) 
—_ . — (1.3.1) 
m+r+1 m 
(m =0,1,...; 7 =0, 4,...). 


1.4 Interpolation durch rationale Funktionen 
Das Problem, die rationalen Funktionen 


2 





p= 
a An, p xP 
o,1,...,2n—1,a(*) = = (1.4.1) 
prgine 
fiir welche 
Po,1,...,2n—1,0(%) =f, (r =0,1,...,2" — 1, a) (1.4.2) 
ist, und 
p=n+1 : p 
Dan, p * 
Bo,1,...,2n,6(%) == (1.4.3) 
An, p xP 
fiir welche p=o * 
Do,1,..., 2n,b(%-) =f, (7 =0,1,..., 2m, b) (1.4.4) 


ist, zu konstruieren, wird durch Anwendung der Thieleschen Interpolations- 
Formeln [5], [6] gelést. Die Folgen reziproker Differenzen 0,,,,(%o, %1, +++» Xn» %q)> 
die man in dem Schema 


Qo (Xo) 
9 (%1) 


2% (*2) 


01 (Xo, *;) 


O2(Xp, %1, Xe) 
01 (Xo, X2) ° 
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anordnen kann, werden nach den Beziehungen 
01(%,) =O Qo(%,)=f/, (7 =0,1,..-) (1.4.5) 


0, 41(Xo, *, AAG x,y, %q) er 0,—-1(%o, 4, adil %,—1) + 


aa (1.4.6) 


Or (Xo, %15 +++) ¥p—a» %a) — Or (Xs %45 -- +s ¥p—gs5 ¥y) 


(ry =0,1,...; @a=r+i1,7+2...) 





gebildet. Die gewiinschten Funktionen 9 1... ¢n-1,a(*),%o,1,...,2n,a(%) sind bzw. . 
die n-ten und (+ 1)-ten Naherungsbriiche des Kettenbruches 


% — Xo %— %, ie 
O1(%, %)+ @2(%o, %1, Xa) — Oo(%o) + 
- Mea ad 4 eeee 
om Orta (For Xp + +++ Xp» Xa) — Op—a (Mos Mp5 «0+ Mp—ay Meta + 





Qo (%) + 
(1.4.7) 





Man errechnet Polynome 4 y,..,,4(*), Bo1,....7,¢(*) durch Anwendung der 
Rekursionsformeln 


Ag 1....,7+1,0(*) = {0,+41(Xo, yy sees My, %q) 7 0,-1(Xo, My eee,) X%y_1, %,41)} x 


(1.4.8) 
X Ag 1....,7,0(*) + (x cS X,) Ao 1....,r-1,a(*) , 


Bo,1,...,7+1,a(*) = {0,41(%o, Hy, oy Xp, %q) — Opa (Ko, Xy,--+, %p—1, %r41)} x 


x Bo 1,...,r,a(%) + (*% — %,) Bos... --1,0(*) (1.4.9) 


aus den Anfangswerten 


A4(%)=1, Ao(x) =f; B_(x)=0, Bo(x)=1 (1.4.10) 
und danach 








Ay 1,...,2%—-1,@ 
Po, 1,...,2n—1,0(*) os z. Stal , (1.4.11) 
A n,a 
o,1,...,2n,a(%) = pe ; (4.4.42) 


Die ,,zentralen‘’ rationalen Funktionen 41... m+2n(*) Om, m+1,..,m+2an41(%)> 
die dieselbe Form wie (1.4.1) bzw. (1.4.3) haben und fiir die 


Pn, m+1,..., m+2n(*7) =f, (r= m,m+1,...,m+2n), (1.4.13) 
Pm, m+, ..., m+2n-+1 (Xp) =f, (7 =m,m-+1,...,m+ 2n +1) (1.4.14) 


ist, sind auf eine 4hnliche Weise konstruiert. Man bildet die reziproken Differenzen 
O,(Xms Xm4is> +++) Xm4,) Mittels der Beziehungen 


Q-1(%) =O Qo(%,) =f, 


Ora (Xmes Xentas ++» Xetra) = Op—a(Xmtar Mme -++> Xmty) + (1.4.15) 


¥m+r+1— *m 
Or (%mtis ¥mtar-++> ¥mteti) — Or(%ms Xmtiy <--> m4 or) 
Numer. Math. Bd. 2 12 
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und danach die u-ten bzw. (n+ 1)-ten Naherungsbriiche des Kettenbruches 








%— X %— Hy 7 
Go (0) + 01(%, %1) + Oe (%o, %1» ¥2) — Qo(%) + (1.4.16) 
ee 4 — Hy eee 
Orta(%os %ys «++» Ft eta) — Or—a(%1) %as +++ He) + 


2.0. Rationale Ausdriicke 

Die erste Untersuchung dieser Arbeit besteht darin, einen Algorithmus zu 
finden, der die GréBen 1... en-1,a(*), 9o,1,...,2n,a(*) in Beziehung zueinander 
setzt und so deren Berechnung erméglicht ohne die Notwendigkeit, reziproke 
Differenzen mittels der Gleichungen (1.4.6) zu bilden, und ohne die Zwillings- 
Rekursionen (1.4.8) und (1.4.9) und die nachfolgenden Divisionen (1.4.11), (4.4.12) 
in Anspruch zu nehmen. 

Bevor dies ausgefiihrt werden kann, miissen explizite Ausdriicke fiir 
Po, 1,...,2n—1,a(%) und J 1. 2n,a(*) abgeleitet werden. Man betrachte den Quo- 
tienten 

Ifo, % — %, fa(¥% — %q), (¥—%5)?, fg(¥% — %4)*, -.-, (% — %an-1)", fal(¥ — %o)"| . (2.0.1) 

|19, ¥ — %, fg(¥ — %2), (%¥— 4 )*, fa(*% — %4)*, .-., (% — ¥an—1)™, fa(¥ — %a)"| 
Die Determinanten in diesem Ausdruck und soweit mdéglich in weiteren Aus- 
driicken dieser Arbeit werden durch die Elemente ihrer Hauptdiagonalen ange- 
deutet. Die Form jedes Elements deutet die Form der Elemente in Spalte an, 
die Indizes jedes Elements deuten die Indizes der Elemente in der entsprechenden 
Zeile an. Die GréBe 1, in der Determinante im Nenner des obenstehenden Aus- 
drucks impliziert, daB jedes Element in der ersten Spalte Eins ist (und daB jedes 
Element in der ersten Zeile den Index Null hat). 

Setzt man x = x, und entwickelt man jede Determinante nach den Elementen 
und algebraischen Komplementen der ersten Reihe, so sieht man, daB der Aus- 
druck (2.0.1) die Werte /, fiir x =x,, r=0,1,...,2”—1,a annimmt. 

Wenn man weiterhin die Zahler und Nenner mit der Determinante 





1 0 0 0 0 0 0 
0 1 0 0 0 0 0 
—x 0O 1 0 0 0 0 
Oo —-x O 1 0 0 0 
2 0 —2x 0 1 0 0 
0 x* 0 —2x O 1 0 
—# 0 3x7 0 —3x O 1 








multipliziert, sieht man, daB der Ausdruck (2.0.1) zu 
fo» % — %, fa %g, %3(% — %), fy ¥4,--., M8n21(¥% — %on-1), fa %2| (2.0 2) 
|19,¥ — %, fe %a, %3(% — Xs), fg X49, -.., Mn ta (% — Xen): fa *2| 
fihrt. Der Zahler von (2.0.2) ist offensichtlich die Summe der Produkte von n 
nicht mehr weiter reduzierbaren linearen Funktionen von x, und ist so ein 
Polynom n-ten Grades in x, wie auch der Nenner es ist. Somit ist der Ausdruck 
(2.0.1) mit der Funktion 1. ¢n—1,4(%) identisch. 
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Auf ahnliche Weise kann gezeigt werden, daB 


9o,1,..., an, a(*) 


= \fo, ¥ — %, fa(% — %a), (¥ — %y)*, fa (*% — %)?, ..-s fan (*% — Fan)”, (%¥ — %q)*t | (2.0.3) 
|19, 4 — 4%, fg(¥ — %g), (% — %)?, fg (* — %)*,..-, fan(¥ — Fon)”, (~ — %4)"+3| 





ist. 
3.0. Drei Hilfssatze aus der Determinantentheorie 


Um allgemeinere Beziehungen zwischen den Funktionen 9 ; ... 2n—1,4(%) und 
Bo 1,...,2na(*) Zu erhalten, wird auf drei Resultate der Determinantentheorie 
zuriickgegriffen, die hier der Vollstandigkeit halber angefiihrt werden. Die beiden 
ersten sind, nach SCHWEINS [7] 


|@y, bg, Cy, «++» Yama» Zn] _ ae |a,, hp| dy |@, he, bg| |d1, Cal ea 


|fty, Dg, Cy, --+> Vn—1» 2n| ~ hy |hy, bg| |hy, bg| | hy, bg, g| 




















(3.0.1) 
ra 1Ba, Ra, Ons >>>» Val | Ors Coe >>> > Ya—e» Sa—al 
4 [hy Da, Cy, -++ + Vn—1| |My» De, Cg, --+ > Zn! 
un 
|a,, De, Ca, +++» Vas» Sn) =a — a, b, = |ae, bg| |b,, Co| les ie 
lbs, Cys +++» In—1» Sah bs bz |be, ¢| (3.0.2) 


___|@g, bg, .--. Yn| [Oy Cg, -+-» Yn—g» 2n—y| 
|bg, Cg, ---, Yn—a/ [Da, Cg,» Vn—a» Zn 


und die dritte [8], die man aus einer extensionalen Identitat, beruhend auf 











|, b2| =| |al, |b21| (3.0.3) 
ableiten kann, lautet 
|a,, be, soe Snr Yntil |r, de, «+1 Zn Mata! 
7? ae a a OF} ey yer (3.0.4) 


oan |a,, be, w+» Sg | |@,, Og, ---> Su» Mnty Sntal 
[ay Dg, --- > Snr Suta| [,, Og, +++» Sus Sutal 


Es sei bemerkt, daB die Resultate (3.0.1) und (3.0.2) in den Formen 

















|a,, bs, .--> Yn—a» 2n| lay, Oar >>+> Ya—al — ba» Mas Bas ++> Hal itis >>> Mamariaaal (3.0.5) 

[Reser esse PamgeSn| Age Ogrsr0s Pama) ae Mer+-+> Fo—al par Oar+++0 Fa—geMel 

| @y,Da,-+-»¥n—1»2nl |r, 0g,---) Yn—al _ ___ Ga bg, ++» Yn [Dr -+ +» Yn—2» 2n—1 | (3.0.6) 
|bg, .--) Yn—1» Zn |ba, «+1 Yn—1| |bo, +++ Yn—1» 2nl [Dg, «++» Yaa 


geschrieben werden k6énnen. 


4.0. Ein Algorithmus 


Um die Schreibweise der folgenden Gleichungen zu vereinfachen, werden die 
Abkiirzungen x — x,=a,, y=0,1,... gebraucht. Aus Gleichung (3.0.5) folgt 


0,1,..., 2n+1,s(*) _ Do... 2n,s (x) 








i [for +++» fon On, O87, fanta Wnt! Ifo» M1» +--» fan Mn, ast? | 
a 119, Oy, +--+ fom Om, OF, fonts OB tal 119,0%,---+ fen On, agt?| (4.0.1) 
aad 10, fis Oy, ++ Om fs OF, Oa trl |e, fy, ---, O8n feast? | 

119, 1, ---+ fon OBn, OF, fonts Mtatal [10,1 ---» fon On, agt?| 


12* 
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und aus Gleichung (3.0.4) 


Bo... 2n,s(%) - Bo... an, 2n+1(%) 

fas Joo, fy ty, ---» Wan fs] pus loo, f,%,---» Ons fan+al 
Yat fe Oy ORE, 1g] [tor free. ORT, tong] (4.0.2) 
lao, fp Oy, ---» ORns famta» ts] 10%, fr Or, ---» On| 


logs fa Oye, Ope, 1g| | or fr, ---. Oe, Iental 








Im Falle des e-Algorithmus [9] und beim determinanten-theoretischen Beweis 
des g-d-Algorithmus [10] erzeugt die Anwendung des Resultats (3.0.4) eine 
Determinante héherer Stufe von gleicher Art wie die in den Determinantenaus- 
driicken fiir die betreffende GréBe schon bestehenden Determinanten. Es ist 
jedoch aus Gleichung (4.0.2) klar, daB die erzeugte Determinante héherer Stufe 
in diesem Falle, namlich 


| 1osh, Xe, [3%, oney On» fen+1%n41> ant} (4.0.3) 
mit den Determinantenausdriicken (2.0.1) und (2.0.3) nichts gemein hat. 
Dagegen fiihrt Formel (3.0.6) zum Ziel. 


Es sei definiert 


\fo. O%,  ) agthy, fs agtt| 























ba tered (4.0.4) 
% m+1 
Mot alan in 1 aii 9o,1,..., 2n,s (*) , (4.0.5) 
fe, %,---, agtt,, f, apt? 
| 19 = " mI | = %o,1,...,2n+1,5(%)» (4.0.6) 
» Gay -++5 Me 
le, Gs, ++», fan Gan a t* 
tin a, Ashe | = %,1,...,2n,s(*)» (4.0.7) 
1 1 re _f, ak 
7 —- a 4 =Wo1,...,2n+1,s(%)> (4.0.8) 
los tis --»: fan Ghat, af 
: ev h> ahs a L = Wo,1,...,2n,s(*) « (4.0.9) 
Aus der Gleichung (3.0.6) folgt 
[fs ob, aS TY, «+5 Gans fants a [fs aft, aft}, ..., Xen! 
|agt*, ++) ams fantal Jagr, ..., fen—s Cen—1» %en| 


(4.0.10) 


- Jag, fo OO, +++) %en—a, fen| Ifo apt}, ORT, 2005 Oe n+; | 
Jagt?, ..., oem, fan+al Jag +?, ++) fan—1 Sen—1> Cen| 


Multipliziert man Zahler und Nenner der rechten Seite der Gleichung (4.0.10) 
mit der Determinante 





[for Or» +++s fen AZ| (4.0.11) 
so kann man sehen, daB Gleichung (4.0.10) als 


Oo ,1,...,2n+1,5 (*) = Xo... 2n,s(*) 
(4.0.12) 


Oo,1,..., 2", (¥) Wo, Lease, 2%, 20-+1(%) Sgn4y 
Go, 1,...,2, 2n+1(4) 





geschrieben werden kann. 
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Aus der Gleichung (3.0.6) folgt 


Ifs ag +}, apt}, e+) hon entail ou) Ifs ag tt, apt, AE Oe n | 
logt*, ..., den, lentil |agt?, ..., fen—1 Cen—1» Cen! (4.0.13) 
_ last, fo WG, ---» tan—1» 1am | fo %o*?, aft, ..., Genta 
JaStt, ..., fen—s Man—1» Xan| Jagtt, ..., den, Iental 








Multipliziert man Zahler und Nenner der rechten Seite der Gleichung (4.0.13) 


mit der Determinante 
| 4qnttgs>>s Fen Ment (4.0.14) 


so sieht man, daB Gleichung (4.0.13) als 


Mo,1,..., 2n-+1,s (*) par XsWo1,..., 20,5 \%) 
Mo, 1,..., 2,5 (¥) o,1,...,2", 2n-+1(%) Senta (4.0.15) 
Mo, 1,..., 2”, 2n-+1(*) 





geschrieben werden kann. 
Aus der Gleichung (3.0.6) folgt 


Ifs as, 26, fp at *, ..., tan—a, fen—a> ton| fs ae, ao, fy n*, 4, en—2,fen—r| 
las, A af %, ..., an—as fem—a» ten| [a0 frat *s «+ Cam—a» fan—a| (4.0.16) 
las, foao—*,.--.fan—a» 1en—11 Ifo OO, OT, .--, Xen—a, fen| 
[og fat", ...,@en—o fem—a» fonl 100, fat 7, -.-, ten—a, fen—al 








Multipliziert man Zahler und Nenner der rechten Seite der Gleichung (4.0.16) 
mit der Determinante 


| 49, far Qgs+>>sfan—1%en—al (4.0.17) 


so kann man iiberdies sehen, daB sich Gleichung (4.0.16) als 





o,1,...,2"—1,8(*) 9,1, ...,2n—1, an (*) 
Wo,1,...,2n,s(*) — %,1,...,2n—1,s(*%) = — “4 = — =) A (4.0.18) 


schreiben 1aBt. 
Dividiert man schlieBlich Gleichung (4.0.2) durch Gleichung (4.0.1) so erhalt 
man 


Bo,1,...,2,s(*%) — 9o,1,..., an, 2nt1 (4) 
Po,1,...,2n+1,5(*) — O94, ..., an,s (*) 














4.0.19 
Hor tr, ++ Re tI fs OFT [Oo fro, +» oe *| — __ %o1,...,ant+1,5(*) ( ) 
|19, Oy, se asd te| |, fy %,--., 08d, fs ag tt| Xs o,1,..., 20, ¢(*) 
Auf eine ahnliche Weise kann gezeigt werden, daB 
Oo,1,..., 2n,s(*) =—s %;Wo1,..., 2n—1,5 (*) 
(4.0.20) 
_ Oo,1,...,2n—1,5(*) ®o,1,...,2n—1,2"(%) Sen 
Oo,1,...,2n—1, 20 (*) , 
7 ,1,. .,2n,s(%) —s X>Wo1,..., 2n—1,s(%) 
(4.0.21) 


_ %o,1,...,2%—1,5(*) o,1,..., 20-1, 2"(%) Sen 





%o,1,...,2"—1, 20 (*) , 
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Wo 1,..., 2n-+1,: (*) aa Oo,1,...,2n,s(*) 





4.0.22 
— — St wan, $(¥) 9,1, ..., 2m, 2n+1(*) ( 
o,1, ..., 2", 2n-+1(*) ’ 
o,1,...,2n—1,8(%) — Po,1,..,2n—2,n(%) ay, vey 20,8 (*) (4.0.23) 





Bo, 1,..., 20, s(¥) — Po,1,..., 2n—1,s (*) Os Wo, 1,..., 2n—1,8(*) 


4.1 Vereinfachung der Formeln 


Nun gilt jedoch, daB, obwohl die Formen der Determinantenausdriicke fiir 
Funktionen mit gerader und mit ungerader Anzahl der Argumente verschieden 
sind, deren Ausdriicke denselben Beziehungen geniigen. Schreibt man nadmlich 


0o,1,...,2n+1,5(*) =O8h41,  90,1,..., 20,2 (*%) = Ahn, (4.4.1) 
Mo,1,...,2n+1,s(*) = En 41 Mo,1,..., 2n,s(%) = m$),, (4.1.2) 
o,1,...,2n+1,s(%) = oft, (o,1,...,2n,s(%) = w§),, (4.1.3) 
Po,1,...,20-1,8(*%) =*h42, 0,1, an,5(%) = Mh, - (4.1.4) 


so werden die Gleichungen (4.0.12), (4.0.15), (4.0.18), (4.0.19), (4.0.20), (4.0.21), 
(4.0.22) und (4.0.23) zu 


— 
(s) of), wf 4, 








ae 
0; = 4,0,21 — oti . (4.1.5) 
r—1 
(s) (r) 
OF oe (s) Ty— 4 Wy—1 Sy 
I, bs a, Wy-1 gee r on , (4.1 6) 
My-1 
(7) (s) 
(s) __ (s) O7—1 Wr-1 
0, = 6,1 — —— (4.1.7) 
Wy-1 
a ,—, aff) (4.1.8) 
9¢'S) — g¢(S) ~~ a. qs) * ane 
r r—1 sr 


Aus den Definitionen (2.0.1), (2.0.3), (4.0.4), (4.0.5), (4.0.6) und (4.0.7) folgt, daB 


(s) (7 
o} xy 4 


OOF coe 
c'= 
r ms) 





(4.1.9) 


4.2 Ein Schema fiir den Algorithmus 


Die Formeln (4.1.5), (4.1.6), (4.1.7) werden nacheinander gebraucht, um die 
GréBen o!), x), wv) zu bestimmen; Formel (4.1.9) wird gebraucht, um die 
rationalen Funktionen x‘) zu berechnen und Formel (4.1.8) wird zur laufenden 
Kontrolle angewandt. In der Tat ist es wirtschaftlicher in bezug auf die Anzahl 
der Rechenschritte, die GréBen o!*), 2‘, («, w{) und x!) als fundamentale GréBen 
zu nehmen. Formeln (4.1.5) bis (4.1.8) werden dann zu 

o) = (a, a4) — OF) 1 (a, wy" 3) (4.2.4) 


Xs Wy—1 a”, ’ 








(s) (r) 

) s mS) y (a, wi”? 1) 

mu) = (a, w'*),) — = — 
My—1 


(4.2.2) 
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$s 
(a, ow) =a, 0), — Oty OY? 1 (Xs ws” 1) 
AX; Wy ) = Hs Oy (7) ’ 
(a, wy” 1) 
WM — 4, a 
= —-. 
ay) — ss (as wfs)) 


Die Anfangswerte fiir das Rekursionssystem (4.1.9) bis (4.2.4) sind 
w=1, (a,m%,)=a,, «% =o%} =/,. 
Die betreffenden GréBen kénnen in dem Schema 


(0) 


at, wo), 
0, 

o), 
oe), 
c) 

a, ow, a, wf)? 
a), ae) 
of!) of?) 
2), xl) 
ct) c(t) 


(s) (s) ee (s) 





a,@'') a, Wy... AW, 
mS), a)... of) 
oa), of)... of) 
x er 
c), eo .. & 

angeordnet werden, wo 
cls) _ (a; ws?) (xf) — sf 1) 
r 


mG — a) 
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(4.2.3) 


(4.2.4) 


(4.2.5) 


(4.2.6) 


ist und den Wert Eins haben sollte fiir alle Werte von s und 7. Unter der An- 
nahme, daB die 4 Argument-Funktionswerte-Paare x,,/, s=0,1,...,4—1 ge- 
geben sind, sind die / Glieder der ersten Spalte des obigen Schemas durch An- 
wendung von Gleichung (4.2.5) berechnet worden; danach wurden die iibrigen 
h—1 Spalten durch Anwendung der Gleichungen (4.1.9) bis (4.2.4) fiir s=7r-+41, 


r+2,...,h—1;7=0,1,...,4—1 berechnet. 


4.3 Die Formulierung in ALGOL 


Die angegebenen numerischen Ergebnisse sind mit Hilfe des provisorischen 
ALGOL Ubersetzers des Z-22 in Mainz berechnet worden. Das komplette 


ALGOL [11] Programm fiir den Algorithmus wiirde lauten: 
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array (alphaomega, pi, sigma, kappa, Cp _1:4-1,,-2, alphag.,_1); 
for s:=0(1)4—1; 
begin pi, _1: =1; alphaomega, _,: = x arg — x,; 
kappa, _1: =/,; sigma, _,: =/,; alpha,: = x arg — x,; end; 
for vr: =O0(1)h—1; 
begin a: = alphaomega, ,_,/sigma, ,_;; 
b: = alphaomega, ,_,/Pi, ,1; 
d: =alpha,/a; e: = kappa, ,_,; 
begin for s: =r+1(1)h—1; 
alphaomega, ,: = alpha, x sigma, ,_, — d x alphaomega, ,_,; 
pi, ,: = alphaomega, ,_, — 0 X pi, ,-1; 
sigma, ,: = alphaomega, ,_, — @ X sigma, ,_1; 
kappa, ,: = sigma, 7-1 % e/Pi, y—1; 
Coy = alphaomega, , x (kappa, ,_1 mae e)/(Pi,,, x (kappa, , _ kappa, r-1)) , 
output alphaomega, ,, pi, ,, Sigma, ,, kappa, ,, ¢, ,; 
end; end; stop; 


Bemerkungen 


1. Es ist angenommen, daB die GréBen h, x arg, und die Vektoren x,,/, . 
s=0,1,...,4—41 schon im Speicher stehen. 


2. Da in den ,,array declarations‘‘ nur rechteckige arrays gestattet sind und 
das angegebene Schema ein dreieckiger array ist, impliziert die gegebene ,,array 
declaration*‘ eine Speicherverschwendung. Notfalls kann man das durch An- 
wendung eines Abbildungsfunktionunterprogramms, das einen dreieckigen array 
bildet, vermeiden. 


3. Uber ,output‘’ Prozeduren ist bis jetzt noch nichts endgiiltig vereinbart 
worden. Die spezielle Form des ,,output‘‘ in dem angegebenen Programm hangt 
daher von der angewandten Rechenanlage ab. 


4. Die Form des ,,output‘‘-Befehls hangt natiirlich von den Umstanden ab, 
unter denen die Interpolation durchgefiihrt worden ist. In einer provisorischen 
Untersuchung wiirde man die GréBen, die im ganzen Schema stehen, ausdrucken. 
In einem schon griindlich untersuchten Fall wiirde man nur die GréBe kappa,_, ,—2 
ausdrucken. 


4.4 Vergleich mit der Thieleschen Interpolation 


Die Berechnung des Faktors (a, «!"),)/o!,, («,a!;)/a,, «, o!7;/(a, w!”;) 
verlangt insgesamt drei Divisionen pro Spalte im obigen Schema. Danach ver- 
langt die Berechnung der GréBen o!), 2°), x!5) und (a, @!*) durch die Gleichungen 
(4.1.9) bis (4.2.3) insgesamt fiinf Multiplikationen, drei Subtraktionen und eine 
Division. 

Die Berechnung der reziproken Differenz 0,,,(%9, %1,..-, %,, %;) durch Glei- 
chung (1.4.6) verlangt drei Subtraktionen und eine Division; die Berechnung der 
GréBen Ag y .. ¢41,s(%), Bo 1... ¢41,s(¥) durch die Zwillingsrekursionsformeln (1.4.8) 
und (1.4.9) verlangt vier Multiplikationen und zwei Subtraktionen, schlieBlich 
fordert die Berechnung der GréBe 1, on—1,5(*) DZW. Do 1... on,s(%) eine weitere 
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Division, d.h. insgesamt vier Multiplikationen, fiinf Subtraktionen und zwei 
Divisionen. Also ist ersichtlich, daB der Gebrauch der algorithmischen Glei- 
chungen (4.0.32) bis (4.0.35) hinsichtlich der Anzahl der Rechenschritte gegeniiber 
der tiblichen Anwendung der Thieleschen Interpolationsformeln vorteilhafter ist. 
Ein weiterer Vorteil ist, daB man die Gleichung (4.2.4) als unabhangige Kontrolle 
nach Wunsch gebrauchen kann. 


4.5 Kontrollen und Invarianten 


Eine weitere Kontrolle gewinnt man durch Anwendung der Gleichung (3.0.2). 
Zum Beispiel folgt aus Gleichung (3.0.2) 


l/s af, @G, fp at", ..-,%en—a» fan—a» 1en| ~ 














}ag, f, On 7, ..., Ogn—a, fen-1> Ian 
at lfs ag, oS, fpaf—*,...,0bn—a.fen—s%an—3, %an—al 
|06, f, Og *, ..., Ohn—e» fen—s San—s» Sgn—sl (4.5.1) 
i oe la, fo XO *, -.-, tan—1» fan—al [fo AO» Ot» ---» Man—a| _ 
|ag, h a, >009 Gign—gs Jan—a) 1a h af—?, oso» Ggn—gl 
— | ag, fo ag}, -++>fan—a> 1en—1| | fo OG, OF, ..-, Sen—1> fen| 
| af, h “a, Sy Sr 1en| | oft, hy - conyp Lon—e@> Ien—1| 


Wenn man Zahler und Nenner des ersten Gliedes der rechten Seite der Gleichung 
(4.5.1) mit der Determinante 


| fora» ++-rfan—2%en—sl (4.5.2) 


sowie Zahler und Nenner des zweiten Gliedes der rechten Seite der Gleichung 
(4.5.1) mit der Determinante 


[49,hrs+++>fan—1 Ona] (4.5.3) 


multipliziert, so kann man sehen, daB Gleichung (4.5.1) als 


Wo,1,...,2n,s(%) — &,Wo 1... 2n—2,5 





- — ei ®o,1, ...,2n—2,2n—1(%) %,1,..., 2n—2,8 (¥) Ps (4.5.4) 
9o,1,...,2n—2, 2n—1(%) 


50,1, ..., 2"—1, 90 (¥) Wo, 1, ..., 2n—1, 5 (¥) 





o,1,...,2n—1, 2n(¥) 
geschrieben werden kann. 


Aus der Gleichung (3.0.2) folgt 








fs agt?, apt, +++» Hon, 1on+1| os 
Jagtt, ..., ten» tom+al 
id Ifs att? aBs, ..0 fen—s On—3, tn—2, fan—1 Og n—1| 
| apt?, --+ofan—s On—s, na, fen—a Xen—y| (4.5.5) 





— latt?, fo OG, ---»fan—s %an—2» ten—al |fo eet, aft’, ...) fan en ve 
logt?, ...,fen—1%en—a» Caml |2OT?, .--s fen—a %en—al 
ae latt?, foae,---»fan—2%an—2 %an—1) 19] [fo%t*, Ot, ---s fen Mans Senta! 
Jogt?, ..., fem—1%en—a» Sans tential lOot?, ..., fen—2 ten—1» Con| 
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Wenn man Ziahler und Nenner des ersten Gliedes der rechten Seite von Glei- 
chung (4.5.5) mit der Determinante 


| Ser tar>s+>han—aMten—sl (4.5.6) 


und Zahler und Nenner des zweiten Gliedes der rechten Seite von Gleichung 
(4.5.5) mit der Determinante 


[19,01,--»,fen®en| (4.5.7) 


multipliziert, so sieht man, daB Gleichung (4.5.5) die Form 


Mo,1,...,2n+1(*) — 4,99 1...., 2n—1,s (*) 








a o,1,...,2n—1,8(¥) %o,1,...,20(¥) Santa Mo, 1,...,2n,s(¥) Mo,2,..., 2, an+1(*) (4.5.8) 
po, 1,...,2"—1, 2n (*) Mo, 1,..., 2, 2n+1(*) 
annimmt. 
Auf eine ahnliche Weise gewinnt man 
Oo,1,..., 2n-+1,s (*) = 0%;09,1,...,2n—1,5(%) 
ae Os Wo, 1,...,2n—1,5(¥) %,1,...,20(*) _ %en+1 9,1, ...,2,s (*) ®o,1,...,an, 2n+1(*) (4.5.9) 
po, 1,...,2n—1,20(*) 90,1, .., 20, 2n+1(*) 


Wieder ergibt sich, daB die Gleichungen (4.5.4), (4.5.8) und (4.5.9) dieselbe Form 
haben, unabhangig davon, ob die Ordnung der betreffenden GréBen gerade oder 
ungerade ist. Man gewinnt schlieBlich 














(7—1) ois) (s) of”) 

18) ses (s) y—y Wy—g Oy-g — Wy—1 O71 
wi) = a,wss) 4.5.10 

? s\7—2 off» wo”, ’ ( 5 ) 
r+1 s“r—1 wo”, arth) ’ 

MH awodt, «SB Of, 092, dy 4 0, wt) (4.5.12) 
rt+1i-” “s“r—-1 ow” oft) 
r—1 r 


Der Gebrauch der vollstandigen Reihe auf der rechten Seite von Gleichung (3.0.2) 
fiihrt zu dem Resultat 


(s) — ls) a? — gl) ar — 
Wo, = OL %, OW} %, 
- - ( 
__ 9 wo?) ay as? ot'tio i wyi—a) of 2 _ @e}-1 Os (4.5.13) 


2r— 27 ’ 
op oR” OP 








s) — -(s) ,~7+1 _ (0) ,7+1 __ 
mh) ss ais O-\%, O_1%, 





S) (1 s (s (27 {s) ,,(27+1 +). 
of ww aap a, W3)_1 O94"? 1 Sart, WE) wrth (4.5.14) 
i @r arti , 

of) w2 1 ee 


s — als) ~ft+i _ (0) ,7+1 __ 
of) =o), af O_%; 





(1 (2 (s) ,,.¢ > 
_ AP aMaeg Os WEP _-1 O21 a 4 Ony WET” (4.5.15) 
2 > 
ow on, of; ¥1) 
a A OD, agate 
2r+1— 4s ee 
(1) »,(S) ~F iy (2 7) (s) (s) ,(27+1) (4.5.16) 
Go Wo %s _ Hay Wey-1 O97-1 Dar O27 





es os... = a 
ow of ayy) 
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+1 __ a9) XH agtt = 


-- (4.5.17) 


oy) w§ af es os w)_ s of7-» = Og y n)_y w2” , 
wow?) oer — mr) ’ 


/ 
mf) = af 








+1 OS); Og at? mol 


of’) = _ a” 





a, 
(1) 6 (7) at (s) (87-1) (s) en. (4.5.18) 
_ Fo Mo es hs WE7—2 OQy—2" Ny OQ7—-1 Wa r— 
wo ofr %— : or) ny 


Bleibt zu bemerken, daB w), bzw. w{), , identisch Null ist, wenn /(x) durch einen 
Ausdruck der Form 


n n+1 

>, 4 #* > a; x8 
= bzw. ee 
dD 5, #8 D 5; #8 
s=0 s=0 


gegeben ist. In diesem Fall sind die Reihen (4.5.13) bis (4.5.15) bzw. (4.5.16) 
bis (4.5.18) Invarianten. 


4.6 Numerische Beispiele 

Dae > 
nd iy! . Die 
Argumente sind durch x,= —2+s, s=0,1,..., 4 gegeben. Als Variable wird x 
x =%3 genommen. In der Tat ist jede GréBe unabhangig, mittels der betreffenden 
Determinationsausdriicke, berechnet worden. Die Resultate sind in Tafel I 
aufgefiihrt. 

Das zweite Beispiel betrifft die Berechnung des Tangens (x), wenn x = 1,5685 ; 
Tangens (x) sei fiir x = 1,566, 1,567, 1,568, 1,569, 1,570 gegeben. Die Resultate 
sind in Tafel II enthalten. Erstens ist zu bemerken, daB die gegebenen Argu- 
mente als x =—4, —2, 0, 2, 4 und der Variablenwert x zu x =1 genommen 
sind. Die Rechnungen sind mit fester Genauigkeit (dreiBig binare Stellen) und 
Gleitkomma-Arithmetik durchgefiihrt worden, und durch die angedeutete lineare 
Argumenttransformation ist das Maximum an méglicher Genauigkeit in der 
nachherigen Rechnung gesichert. In der Tat fiihrt der Gebrauch der zuerst 
gegebenen Argumente zu einem gerechneten Werte von 435,478226 fiir den 
Tangens (1,5685). 

Die Funktionswerte sind einer sieben (dezimal)-stelligen Tafel des Tangens [12] 
entnommen, die angibt, daB Tangens (1,5685) =435,47730. Der mit den Glei- 
chungen (4.1.9) bis (4.2.4) berechnete Wert (435,47732) stimmt mit dem er- 
haltenen Wert durch Anwendung der Thieleschen Interpolationsformeln [13] 
iiberein. Zweitens ist zu bemerken, daB die GréBe cS) eine numerische Eins, 
nicht eine absolute Eins, ist, d.h. es tritt ein Genauigkeitsverlust an bei der 
Differenzen (x{5),—x!),) bzw. (x‘s)—x!S),) durch Ausléschung ein. Ein weiterer 
Genauigkeitsverlust deutet Instabilitat im Algorithmus an. 

Das dritte Beispiel betrifft die Berechnung der n-ten Nullstellen 7),, der 
Besselschen Funktion J,(x) wenn die ersten Nullstellen gegeben sind. Tafel III, 
die die ersten acht Nullstellen mit Differenzen enthalt, zeigt ganz klar, daB 
Extrapolation mittels eines Polynoms in  ungeeignet ist. 


Das erste Beispiel betrifft die rationale Funktion /(x) = 





5,0 

1,0 
208,491 28 
208,491 28 

3,0 

1,0 
263,411 25 
263,411 25 


a, =3 


1,0 

1,0 
357,611 06 
357,611 06 


—1,0 

1,0 
556,690 98 
556,690 98 


— 3,0 

1,0 
1255,765 59 
1255,765 59 


a=—3 
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Tafel I 
a, wy = 1 
a= —1 
oy) = — 5/3 
x = 3 
ob =1 
ty Wy) = —16/5 a,wf = —24 
n= —2 a? — — 26/5 
of = 2/3 of = — 14/5 
(2) — — 3/5 xP? = 21/13 
ci) — 4 c= 4 
4 Wy = — 27/5 tg WE’ = — 36 tg wg = 6 
ng’ =—3 ny’ = — 42/5 mg’ = 36/13 
of) = o=—18/5 of = —36/7 
ngr=—9/5 x =9/7 xg) —3 
cP = 4 cP =1 cf = 4 
a, o——16/5 aw%=—24 a wi=32/5 a, wo) = 48/5 
ny =—4 n= —36/5 ay=120/13 ax?=—68/5 
on =4/3 of =—12/5 of'=—24/7 of =—12/5 
x) = — 3/5 x = 4 xv) = — 3/5 x) — 9/17 
c= c= 1 c= 1 co =1 
a,a)—=—7/5 awP=—-12 aowf=4 a, w= 12 a,0f?=0 
mo =—5 nf = —32/5 af=228/13 ay’=—34 ny=—12 
of) = — 1/3 of) = — 8/5 of = 12/7 os =6 o, =—12 
xq = 3/25 x?) = 3/4 xy) = 3/19 xg = 9/17 xq =9/17 
cS) — 4 cS) = 4 P= cP —4 cP =1 
Tafel II 
164,759910 
—2,0 
— 3,31708076 
345,791 205 
0,999 999 994 
149,119778 + 1,43043923 
—4,0 — 180,400042 
—7,57616345 —227,189421 
394,891 001 435,477192 
1,00000003 1,0 
— 348,199698 — 6,68025803 0,002 21633911 
— 6,0 — 842,479430 — 13,360498 3 
—14,3504619 —1060,98934 — 13,3605020 
498,657 694 435,477437 435,477 314 
1,0 0,999 999 997 0,996 125359 
— 3141,82991 — 90,414 5890 — 0,060 5468750 —3,32978630 
— 8,0 — 3800,862 55 —120,552624 — 0,040 548 7064 
—33,1155421 —4786,67696 —120,552658 — 0,040 548 7065 
863,037 716 435,477676 435,477 313 435,477 315 
0,999 999 994 0,999 999 998 1,00001 33 0,746 594440 
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Tafel III 
n Jo,n 
1 2,404 826 
+ 3,115252 
y 5,520078 + 0,018 398 
3,133650 — 14238 
3 8,653 728 0,004 156 + 11659 
3,137 806 2579 — 9894 
4 11,791 584 0,001 577 + 1765 
3,139 383 — 814 
5 14,930918 + 0,000 763 
+ 3,140146 
6 18,571 064 


Tafel IV zeigt die GréBen x‘), die aus den gegebenen Funktionswerten und 
dem Argumentwert » = 11 hergeleitet sind. 


Tafel IV 
wt np op Xp xp 

2,404 826 

5,520078 33,55735 

8,653728  33,64934 34,405 25 
11,791534  33,69385 34,18119 33,77625 
14,930918  33,72006 34,05062 33,77611 33,775 58 
18,071064  33,73730 33,965 19 33,77602 33,335 66 33,775 84 


MmPPWNHRe OS 


Der Wert von j,,, auf fiinf Dezimalstellen abgerundet ist 33,77582. Der 
Wert von j,,,, fiir groBe m ist durch die asymptotische Reihe [14] 


<ee es. | 0,050661 _ 0,053041 0,262051 OARS 
jon = {n 2 a 4n—1 (4n — 1)8 + (4n — 1)5 } (4.6.1) 





gegeben. Das letzte numerische Beispiel betrifft deshalb eine Art der Stieltjes- 
schen Summationsmethode. Eine weitere Bemerkung iiber diese funktionen- 
theoretische Eigenschaft des Algorithmus wird in einem spateren Abschnitt ge- 
geben werden. 


4.7 Singulare Regeln 


Offenbar bricht der Gebrauch der Formeln (4.1.9) bis (4.2.4) zusammen, wenn 
zufallig eine oder mehrere der GréBen o),, 2!),, w!,, x, gleich Null sind. 
Diese Schwierigkeit kann man iiberwinden durch einen Umtausch der Argument- 
anordnung. Diese MaBnahme zerstért jedoch die eventuell vorhandene Eigen- 
schaft des x, einen Grenzwert zu besitzen, es sei denn, daB die Argument- 
anordnung vollstandig umgekehrt wird. Die erwaihnte Schwierigkeit besteht fiir 
fo=x, =o, =0. In diesem Fall ist es einfacher zu jedem Funktionswert eine 
Konstante zu addieren, und nachher dieselbe Konstante von dem berechneten 
Wert von /(x) zu subtrahieren. Als Beispiel betrachte man die Berechnung des 
Funktionswerts von sn(x|0,2) [15] fiir x =0,3, wenn fiir x =0,0; 0,1; 0,2; 0,4; 
0,5; 0,6 Funktionswerte gegeben sind. Hier ist sn (0,0| 0,2) =0 aber aus dem 
angegebenen Grund ist jeder Funktionswert kiinstlich um Eins vergréBert worden. 
Das Resultat sieht man in Tafel V: 
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Tafel V 

wo x? Xf) nf) nf 
1,00000 
1,099 80 1,2994 
1,198 41 1,29761 1,295 85 
1,387 52 1,290 64 1,293 50 1,29497 
1,47595 1,285 57 1,292 32 1,295 14 1,290 32 
1,55912 1,27956 1,291 14 1,295 26 1,294 68 1,294 68 


Das Ergebnis 0,29468 stimmt auf fiinf Stellen. 

Fiir die Genauigkeit der rachherigen Rechnung ist es am besten, wenn man 
als Konstante — $(/max-+/min) (vorausgesetzt, daB dieser Wert ungleich Null ist) 
addiert. fimax bZW. fmin bedeuten dabei den gréBten bzw. den kleinsten der gege- 
benen Funktionswerte. Die Konstante Eins im obigen Beispiel ist aus Bequem- 
lichkeit gewahlt worden. 


MWPWNHNHO F 


4.8 Vergleich mit der Polynom-Interpolation 


Bei den bisherigen Beispielen (Tangens (x), 7o,,) war die Interpolation mit 
Hilfe der Aitkenschen Methode unpassend. Das letzte Beispiel ist jedoch von 
MILNE-THOMPSON [16] gegeben, um die Aitken-Methode zu erlautern. Die ersten 
zwei Spalten des Aitkenschen Schemas sind namlich dieselben wie in Tafel V, 
die iibrigen vier Spalten findet man in Tafel VI: 


Tafel VI 
0,295 83 
0,29356 0,29469 
0,29428 0,29471 0,294 67 
0,29466 0,294 74 0,294 67 0.294 68 


Man bemerkt, daB die numerischen Ergebnisse, die mit Hilfe der Polynom- 
Interpolationsmethode hergeleitet wurden, etwas besser als diejenigen sind, die 
mit der rationalen Funktions-Interpolationsmethode hergeleitet worden sind. 
Dies ist jedoch kaum verbliiffend, da die betreffende Tafel von sn (x| 0,2) fiir 
Polynom-Interpolation bestimmt war. 


4.9 Gerade Anzahl der symmetrisch verteilten Argumente 


Wenn die gegebenen Argumente der Reihe nach —a, a, —b, b, —c, c, ... sind, 
kann man eine wesentliche Verminderung des Rechnungsumfangs durch Anwen- 
dung eines Kunstgriffs von AITKEN [17] erreichen. 


Erstens sei angemerkt, daB 








f(y) *-y 
f(x; 9) = /2y (4.9.1) 
f(—y) x+y 
eine gerade Funktion von y ist, und zweitens, daB 
f(x, x*)=f(x) x*x=—a,a,—),),.... (4.9.2) 
Man nimmt deshalb als Anfangsfunktionswerte /(x; a”), f(x, 6?),..., Anfangs- 


argumente a?, b?..., und als neue Variable x. Auf diese Weise werden zweimal 
so viel Funktionswerte als Anfangswerte beriicksichtigt. insgesamt ist die 
Rechenarbeit bei dem neuen Schema ungefahr ein Viertel von der, die bei 
der normalen Interpolation anfallt. 





oy 





of 


s) 
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Als Beispiel diene die Berechnung der Jacobischen Zetafunktion [18] Z (x/| 0,6) 
fiir x =0,11; die Funktionswerte sind dabei fiir die Argumente x —0,00, 0,04, 
0,08, 0,12, 0,16, 0,20 gegeben. Hier ist /,=0, und deshalb ist wieder eine additive 
Konstante (gleich eins) kiinstlich hinzugefiigt worden. Die Ergebnisse der nor- 
malen Interpolation werden in Tafel VII gezeigt. 


Tafel VII 
xls), xs) x{s) xis) xfs) fs) 
1,0000000 
1,0133469 1,03567040 
1,0266172 1,0365986 1,0365200 
1,0397350 1,0364237 1,0365485 1,0364775 
1,0326262 1,0361805 1,0363968 1,0364875 1,0364739 
1,0652186 1,0358702 1,0363345 1,0364936 1,0364739 1,0364739 


MWUbWNRO 


Durch eine lineare Argumenttransformation bekommt man die Argumente 
—10, —6, —2, +2, +6, +10 und die Variable wird x =1. Nach Anwendung 
der Transformation (4.9.1) und Rechnung in der iiblichen Weise kommt man zu 
den Resultaten in Tafel VIII. 

Tafel VIII 
Sx) xf) x{°) 
0 0,0364 55550 
1 


0,0362 59825 0,0364 73899 
2 0,0358 70230 0,0364 73841 0,0364 73931 


In beiden Fallen ist der berechnete Wert, 0,0364739, fiir Z(0,11|0,6) auf 
sieben Stellen genau. 

Leider ist es nicht mdglich, denselben Kunstgriff anzuwenden, wenn ein 
Wert /, im Zusammenhang mit /,, /_4, fs, f_», ... gegeben ist. Obwohl /(0; 0) =fy 
ist, gilt aber /(x;0) =f)— x fg man kann deshalb nicht die ersten Anfangswerte 
des neuen Schemas berechnen, ohne die Ableitung von /(x) fiir x =O zu kennen. 


4.10 Gebrauch von gegebenen Werten der Ableitungen 
Wenn man in den Determinantenformeln (4.0.4) bis (4.0.9), x,=%+s At 
s=0,1,..., 2” bzw. 2n+41 einsetzt und Az gegen Null strebt, kénnen die gege- 
benen Determinantenausdriicke in folgender Weise umgestaltet werden. 











1 1) 0 ve 0 18) bes 0 1 
XK —1! 0 ees fe) ce) ss’ 0 Os 
ag — 2% 2! yan 0 0 en 0 a? 
ag — 3a§ 6a 0 0 0 as 
aft? —(n+1) af (n+1)maf-?...(—1)"41(m44)! 0 2... 0 att 
’ id ie fmt 1) f+? a” f@ n+ 2) ts 
fa (fay fay”, (fat (Fant... (fa)tMtD fre, 
_ | tet (fam (fa (fame (famintD  (fanfOmt fot 
n+2 f"+ 1) att f@ n+2) 
FS e (fayntn ~~ (fayant™ 
r=0 ; ; 


(fam)n+n Ma (f a%)(2"+2) 


(4.10.1) 
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1 0 18) i) 18) 0 1 
% —1! 0 0 0 0 Os 
ag —2m% 2! 0 0 () a? 
af — 30% 6a 0 0 0 as 
at —nat—? n(m—1)ah—2...(—1)""l 0 0 a” 
f ad wa {™ fern feet” hs 
fa (fa)’ (fa)’” (fa) (fa®)tY ... (fa)e*t® f, a, 
, fam (fam) (fa)... (fam (fa)™D ... (fa OMt™ fat 
of) = Far pant (4.10.2) 
Frn| damtm - ayann 
hi (f at —1)in+ 1) Le (f ot—1)(2 n+1) 
Xo -—7 0 0 0 0 Os 
az — 2% 2! 0 fe) (8) a2 
aft! —(n41)of (n+1)mat-?...(—1)"+4(m +4)! 0 0 aft 
, F PO in ee Oe eS 
fa (fay’ fay” Fayre — (Fay... (Fa)*™ fray 
« —tfam fem)’ fat” (fant faryet... fae)” feat 
Os), rs oy ; . . : (4.10.3) 
af — 2m 2! 0 0 0 
of —nat-) n(n—1)af-? ... (—1)""!l 0 0 
f’ ti fj f*rt 1) {? n) 
fo* (far) = (fa)” (far) (famm+9 . (famyiam 
Xo —1 0 0 0 0 Os 
of —2a% ~ 3! 0 0 0 as 
aft) —(m44)af (m+1) nah? ...(—1)"t1(n 441)! 0 0 ant} 
} f’ ad faty [mr® fieaty hs 
of) a fort (f +2)’ (f t+)” (f gntiyinty jantiynen ™_ (fon+1)an+1) hs atti 
oo+3 Oe —1 0 0 0 0 
att} — (m+ 1) af (m+ 1) mof-) ... (—1)"+1(m +1)! 0 0 
} f’ ye fat f™t+® faery 
fa” — (fam)’ (fa) (fart —(fanynt ... (famyanen 








(4.10.4) 
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1 
Qs 


io) 


18) 
0 


0 
0 


0 
O 


O 
=~ tt 


1 

Qo 
0 ont 

(fa)@*—) f, as 


"41(m41)! 0... 
(faye 


(fayn+0 


oft 


y 


nteonae (n +1) mas-? -(- 


(f a)" (f a)” 








(fat) (fa%)” (fam) (faM)M+2 .. (foM2M—D fat 


fam 





nf, = 


(s) 
T3n+1 


0.3) 





10.4! 


(fa)*+) (f a)"*+2) <= (f «)(2™ 
IT r! ; . ‘ 


r=0 





(f am)(2 n) 


(f any) (f nye c 


oO 
oO 


0 
0 


0 
0 


0 
~¢! 


1 
ry 
0 
(fayem 


0 
(fayint 


n—1 


att —(n +1) a% (n+ 1) nak 
fo (f a)’ (f a)’ 


.(—1)"+2(n 41)! 
(faynen 


fantiynse a fanisyam i. 





fort (f +2)’ (f ot +2)”” (fat +1yin+a) 


(4.10.5) 


1 
as 


attt 


fs ts 


atth 





ion (fayentn l 


IT 1! 


r=0 


| (faynr® (foyer... 


(fam)in +2) ayn ioe fanyant| (4.10.6) 


x{S) = Zahler von (4.6.3)/Zahler von (4.6.5), (4.10.7) 


xs), = Zahler von (4.6.4)/Zahler von (4.6.6). (4.10.8) 


In diesen Ausdriicken ist 


(a= OF 


(4.10.9) 
pao? 


) (— 1)? — me)? #8=2). 


Unter der Annahme, daB nicht nur /,, sondern auch fo, fo’, ..., /?”) gegeben sind, 
kann die (2”-+-1)-te Spalte des Schemas durch unabhingige Berechnung der 
Ausdriicke (4.10.1), (4.10.3), (4.10.5) und (4.10.7) aufgebaut werden, bzw. sind 
fo. fo, ---, {2"*” gegeben, dann kann die (2+ 2)-te Spalte des Schemas durch 
unabhangige Berechnung der Ausdriicke (4.10.2), (4.10.4), (4.10.6) und (4.10.8) 
gebildet werden. 

Auf diese Weise werden nicht nur die Funktionswerte, sondern auch die 
gegebenen Werte von Ableitungen beriicksichtigt. 

Aus Bequemlichkeit sind die Ableitungen in der obigen Darstellung an der 
Stelle x = x, genommen worden, aber es ist klar, daB dieselbe Behandlungsweise 
auch bei allgemeinen Verteilungen angewendet werden kann. 


4.11 Fehler-Fortpflanzung 


Ein Fehler des 6a, in a, bzw. df, in f,, pflanzt sich langs einer Zeile bis 
zur (s+ 1)-ten Spalte fort, und danach iiber den ganzen Teil des Schemas, der 
rechts von der (s+1)-ten Spalte liegt wie in Fig. 1 angedeutet ist. 


Numer. Math. Bd. 2 13 
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t% * Durch partielle Differentiation der Verhalt- 
nisse (4.1.9) bis (4.2.3) nach a, bzw. f, und 
nachherige Multiplikation mit da, bzw. df, 
ty Cue ee kann man ein Fortpflanzungsgesetz des absoluten 





cle Fehlers in den GréBen o), x), (a,c), x!) be- 

kommen. Dieser Kunstgriff ist urspriinglich 

* * * in [10] angegeben worden. 6 {o'); «,} bezeichnet 

den Fehler von o{") wegen des Fehlens da, in «,; 

Pon ee *, entsprechende Ausdriicke bilde man fiir die 

Fig. 1 Fehler in den anderen GréBen; man gewinnt so 
nacheinander : 


db {o!); a} = 6 {(a a, ov! )) a} — ED) 8 f0f2h: a, } 


y¥=0,1,...,8’—1, 


B01; fo} = Af (co; fo} — PEED (0h; te} 
(4.41.2) 
y=0,1,...,8’—14, 


b{x); ; a} = {( ( ow) )); a} oe d {x Ws 5 hy BY 


y¥=0,1,...,8°'—14, 


8 fee"); fo} = 8{(oy- eof) fo} — Cees) 0 fol f.} 
(4.41.4) 
¥=0,1,...,8 —1, 


, ) 
BE of); iy} = Of + ty BOL; 9} — EEA (a Os ay} 
ellis (4.11.5) 
¥=0,14,...,8'—1, 


é ; (), a} 7) Me), - __% O71 5 
{(a,.007°'5 fe} = ay {op fe} (a, wo 3) {(a FO) 1) (4.11.6) 
¥=0,1,...,8°'—4, 


Boel"); 1g} = (oe?) [elds fol); 2} — wl" 5 fae!” 9] 


4.11.7 
y=0,1,...,8 —4, ( 
db {x8*); hy = (x!) -1 x” 0 ao); P — x!) 6 ms). 
a aha 
und nachher 
5 {0}? 5 a. } = (of) [ols 5 {(ce, co!) 4) ; x} + [(a*4) — 0] 6 fo! ;0,-} — — 
11. 

enis 9), d{(a, wo!) 1); ay} — (a, (a, w!”) ») b {ol 1; ay} ], ( 9) 

0 {o!) ;2,.} = (m1!) -1 [a é {(a, w!),); ay! + [ (a, ow), ne!) bin” _ 
ria} + U PU ae 


— mh) b{(a, of); ay} — (a, ws) d {a1 ; ay}, 
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5 {(cr, cn!) ; o1} = (ce, co!” 4) [oe, (ax, c0!"21) 5 {of 15 a5} + [a,o!, — (a, w!*))] x 
x bfx, wo! 15 a} — a, [ol b{(a ol) 1); ay} + (aco!) b{o! 1; a,}]], 


5 foe!) cx} = (at!) 4 [orl 5 fol 5 oc} + 281 b fol”; a} — x) d{a!); ae}] (4.11.12) 


(4.41.14) 


mit ahnlichen Ausdriicken fiir 5 {o!); f,} etc. nur ist in der Formel (4.11.9) bis 
(4.11.42) r=s’, s’+14,...,S=r7+1,7+2,. 
Die Anfangswerte fir das Rekursionssystem (4.11.4) bis (4.11.12) sind 


{0%} ; a} = dfn"); a} = dfx"); a} = 0, (4.14.13) 

J {(x,-co©}) ; x} = day, (4.11.14) 

{x}; f.} = 6{(a-o}); f.} =0, (4.11.45) 
db {0}; f} = be); fy} = Of, (4.41.16) 
fo!) 1; 06} = bf); a} = d foes; ay} = b f(a, ws) ; ae} =O, (4.11.17) 


b fo; fe} = bat; fe} = 0 fe; fe} = O{(a@84) 5 fe} =O, (4.11.18) 


wobei in den Gleichungen (4.11.17) und (4.11.18) s=s’+414, s’+2,... ist. 

Wenn die Fehler 6«,., df, klein sind, bieten die Formeln (4.11.1) bis (4.11.12) 
eine geniigend genaue Beschreibung der Fehlerfortpflanzung. Zwei numerische 
Beispiele folgen, die die Berechnung des Tangens (1,566) wie in Tafel III, be- 
treffen. Im ersten (Tafel VIII) wurde x, als —1,999999 statt —2 genommen, 
im zweiten (Tafel IX) wurde /, als 263,41135 statt 263,4125 gewahlt. 








Tafel VIII 
Y 
4 —1 (8) 1 : = 3 
0 
0 
O 
8) O 
—1,0 —1,36 
18) —1,01 
(@) —1,01 
O 18) +0,0271 
0 0 +2,78 
0 0 +2,41 
2 18) 0 +1,21 
0 0 —0,0638 +17,63208 
(0) (8) +4,18 — 26,00 
(¢) 0 + 4,60 — 23,30 
3 0 0 — 0,61 +0,31 
0 0 —0,577  +301,82 — 3684 600,0 
8) 0 + 5,60 —1,55 +402,949658 
8) Oo +410,50 —1,40 +402,949659 
4 O i@) —1,19 +0,60 +0,22 


13* 
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Tafel IX 
% 
s : —1 fs) 1 2 3 
0 
(0) 
O 
O (0) 
0 +3,01 
Oo 0 
—1,0 —2,40 
1 —1,0 +2,50 
(e) 0 —0,0996 
2 0 0 —2,35 
0 (0) +0,2318 —48,5610962 
(0) 0 — 9,03 +0,849 
¢) 0 —7,90 +0,769 
3g 0 O +1,71 — 0,86 
0 0 —2,093 —874,93896 +853254,3 
0 0 —12,00 +5,10 — 436,767 759 
O (0) — 18,10 +4,63 — 436,767 883 
Tafel X 
Y 
s —1 0 1 2 3 
0 
0 
0 
0 0 
—1,0 +2,36 
O —1,01 
—1,0 —3,41 
1 —1,0 +41,80 
0 (0) —0,0724 
(8) 0 — 3,24 
0 10) — 4,75 
2 0 0 —2,22 
0) rt) +0,1686 —31,8145752 
3 0 —1,12 — 0,54 
0 0 +1,521 —572,81694 +733418,9 
0 0 — 6,40 +3,54 — 285,749357 
4 0 0 +2,24 —1,12 —0,015 





Wenn die Fehler geniigend klein sind, kann man wieder eine obere Schranke 
fiir lie Fehler in den GréBen of), 2°), x‘), (a, w) gewinnen. Anfangswerte fiir 


die obere Grenzen von da), da, .. 


., Of, Of,, ... kennt man. Man wiederholt die 
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obige Analyse mit s’=0, 1, ... Danach ist 


13 fo, zu!" x"), (c,00') ; a9, ty, «+s forfrs +++} 
< Pd fo"); 2, 20", (aco!) soy fe}, (4.41.19) 
s 


wo die Ausdriicke auf der linkeri Seite der Ungleichung (4.11.19) die Fehler in 
ao!) bzw. a) etc. wegen des Totalfehlers da in a, da, in %,..., df) in fy, df, 
in f,,... bedeutet. 

Tafel X zeigt den resultierenden Fehler, wenn man beide friiher erwahnten 
Fehler (Tafel VIII, LX) gleichzeitig begeht. 

In den Tafeln VIII, IX und X sind alle GréBen mit 10° multipliziert worden. 
In der Tat sind die Fehler, die in den drei vorhergehenden Beispielen gemacht 
worden sind, etwas zu groB, weil nur fiir die ersten Spalten die betreffenden 
Gleichungen gelten, da die fortgepflanzten Fehler in der letzten Spalte nicht 
mehr klein sind. Trotzdem sind die drei Beispiele lehrreich insofern, als die 
Berechnung der GréBe x‘*) genau b'eibt, obwohl die Berechnung der GréBen o!, 
x) sehr ungenau ist. In der Tat sind die numerischen Ergebnisse fiir x‘) besser 
wenn beide Fehler gemacht worden sind, aber man sollte dieses Prinzip nicht 
verallgemeinern. 


4.12 Inverse Interpolation 


Das Problem, eine Abschatzung fiir das Argument, das einem gegebenen 
Funktionswert entspricht, zu finden, wird leicht gelést. Man vertauscht einfach 
die Rollen von %), %,..., und fy, f,,.... 


4.13 Numerische Differentiation 


Durch Differentiation der Gleichungssysteme (4.1.5) bis (4.1.9) nach x kann 
man ein Rekursionsschema, das zur Berechnung der Ableitung von /(x) fiihrt, 
gewinnen. Man bemerkt erstens, daB w‘*) unabhangig von x ist. Danach gewinnt 
man 








4 ws) , 
of, — 2 ol, (4.13.4) 
r— 
(sos) WH OF, | WH OF OF") 
O; Wy-1 a”, ao, oi ’ (4.13.2) 
(s)’ (r) _ ss HD, wf, ty coy rps Oy O42 y ya Mya 
N;,’ —W,.1>— im “ae 
1p) = {ah ach + of?) xh, — vel? aff” (oe) 2 (4.13.4) 


und durch Anwendung von (4.13.1) 





, of, wi? a, wh) of, wi?) 
o) —ow",= — . oi, 1 + - . one rd (4.13.5) 
mit 
of), — gpl, — x6 — 9, (4.13.6) 


Ein Strich deutet dabei Differentiation nach x an. 
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Man kann aber dieselbe rationale Approximation der Ableitung von /(*) durch 
Anwendung der Thieleschen reziproken Differenzen gewinnen. Man berechnet 
die Polynome 46 1,....-+1,s(*), Bo,1,...,7+1,s(%) mit Hilfe der Zwillingsrekursions 
formeln 
Ao,1,...,741,8(%) = {Op41(%os %1s «++ Xp» %s) — Opa (Xow ¥as +++ Xp—1» Xp 4a)} X 

Xx Ao... 7,8 (x) + (x ies x,) Ao... »7—1,s (x) + Ao)... r—l1,s (x) , 
Bo... »rt+i,s (x) a {0,41 (Xo, Ay, +++y My, Xs) — Or-1 (%o, Hyp +++, %y_1, %,41)} x 


Xx Bo....,7,s(*) + (x > %,) Bo Rory r—1,s(*) + Bo peasy r—1,s (%) 


in Zusammenhang mit (1.4.8) und (4.4.9) und gewinnt x‘*” aus 
2) “ag (Bos... (x))2 {Aoi oop, @ (x) oe: x\*) Boa Siri 1,8 (x)} ° (4.13.9) 


Dieser letzte ProzeB ist in der Tat hinsichtlich der Zahl der Rechenschritte so 
viel wirtschaftlicher als der Gebrauch der Formeln (4.9.1) bis (4.9.5), daB eine 
numerische Untersuchung der letzteren nicht durchgefiihrt wird. 


(4.13.7) 


(4.13.8) 


5.0. Ein zweiter Algorithmus 


Um eine ,,zentrale‘‘ Art von Algorithmus, analog zu dem Nevilleschen ProzeB 
fiir Interpolationspolynome zu gewinnen, miissen entsprechende GréBen definiert 
werden. Die allgemeinen Betrachtungen wie im ersten Teil des vorherigen 
Kapitels bleiben dieselben. Die fundamentalen GréBen sind also 





| fon» Amtyzs-++s On+2r—1,/m+er On +2 »| __ _(m)’ 
L miarl (my (5.0.1) 
| fn» Amtpssses On+2r—1| 


Lins Synz mtar ar, taiarsal — gl (5.0.2) 


| fm Omtis+s+>tm+er On+2r| 




















es ora om ace wea) 
88 
sores er, ae ane 
dee Ne, = G0 
| fms Smt» +++» Xmt+ar—1» fmtar Omt2rl =o)”, (5.0.7) 


| fmt-1> Star =++> Om+ar| 

mit ahnlicher Definition der GréBen of"), af”, af", wf”, wf").1, (der Strich 

ist jetzt nur ein Unterscheidungszeichen und deutet keine Differentiation an). 
Aus der Gleichung 


lobe ter+a, hm+sar nor +++» mtd, fm| CAs e fmtay Ontar, eee, On+-1| 
lfmtar Oetars +++» Smta> fel lft a9 tars +++» Sets | 








| (5.0.8) 


lfm+arts Ont arti» n+er, tree Om+s>fm+1| |agstite sy, fm+sr—1 On+er—1> ose, Om 
lfm+ar nar see, Amt, hms | \fm+ ar Gnters sees Om +-1 | 
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gewinnt man aber 








off. — Cnt ar+iMbe ilo wee ? (5.0.9) 
und aus der Gleichung 
loc, fot sine» <--> Omi a1, Imi ari —_ loat, fm1 Cntr, +++» omtarl 
lfm-+1 Ota sso Omtsrs imter+al lfm+1 Otis s++> Om+ar| (5.0.40) 


— \fm On » Ont» +++» Om+ear—1> fm+erl Ase hm On » 2295 Oom+ar+a| 
| fa a eoep On+eorlm+ert+il lfm+1 m4» see Om +97] 





bekommt man 



































(m)’? (m+1)° 
” \”’ = Co. Ly @. 
OFF, — yooh” = — =e eta (5.0.14) 
2r 
Auf eine ahnliche Weise gewinnt man 
(m-+1) (m)’? 
, +1)’ _ 6: Lym 
OF a — Om 4274 20Ert a iE a stot ? (5.0.12) 
2r+1 
(m)"? (m+) 
r 7 G 1 %m @ 
of”). — Oy COS, we a See amt Yeti . (5.0.13) 
(m+ 1) (my)’* 
4 sy 1 Ay, 
cL — Oy + 9741 OS =<" af mim)’ sf > (5.0.14) 
27 
” ” my’ (m-+1)y 
) aa 1 om @ 
Te — yy Ey” = — EF atiy Talay (5.0.15) 
(m+ 1) (m)’’ 
, +1) Ki 4 Ly, 
27+1 
” ” (m)’’ (m+1y 
—_ 1 Om @ 
mis — agents = — SEs Omnae re OEY 5.047) 
, (m+1y (my 
+1)’ @ Ad 
OFr 1 — Of) = — Ee (5.0.18) 
2r 
ve - (my’ ~im-+1) 
(m) , Wer 02 
Wor41 sa O02, — “OTF ? (5.0.19) 
» (m 1y (m)’’ 
(m) ie em @ G: 
M2742 — oft? = — - , (5.0.20) 
2r+l 
” ” wm’ om 1y 
OFF -2 — Ofer = — Aart (5.0.24) 
27+1 


Es kommt wieder einmal vor, daB die Form der betreffenden Gleichung (5.0.9) 
bis (5.0.21) nicht von der Anzahl der Argumente abhangt, d.h. daB die Glei- 
chungen (5.0.9) bis (5.0.21) auf folgende Weise gegeben werden kénnen: 











OMY orgy coll? 
ee 
” ve om, 4 ted 
OS" — Oy Oy = — SESH wa , (5.0.23) 
a swe 
ai) — yy 4 vent?) =-—-— i ~ : ) (5.0.24) 





” ed mmr, im- 1 
a — age = — AE cate eL (5.0.25) 
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wi! — gmt’ = — oe, (5.0.26) 
em” — gmt’ — — ee (5.0.27) 
mit den beiden weiteren Gleichungen: 
xt”) — ae ™ eo, (5.0.28) 
oe wom 


Man kann ein Rekursionsschema zwischen den einmal und den zweimal ge- 
strichenen GréBen gewinnen. Zum Beispiel: 








Of” (OL — yg s O84) = OE FOP” + EY (Cy — Ores)» (5-0-30) 
of” — Om4s weet” pt ha 
TY — Chay 4 5 Ey ™ a”, , (5.0.34) 
my’ (mM) 
i) == a (5.0.32) 
i”, _ pA Yat a = Ps lt a (wy = oti") 
~ San * an or = On+s ae (5.0.33) 


aber in jedem Fall ist der Rechnungsumfang dieses zweiten Algorithmus ungefahr 
zweimal so groB wie der des vorigen Algorithmus; daher wird auch hier eine 
numerische Untersuchung nicht durchgefihrt. 


6.0. Konfluente Argumente 


In diesem Abschnitt werden Ausdriicke fiir die GréBen o!, 2‘) etc. hergeleitet, 
wenn die Argumente %, %,,... gegen einen gemeinsamen Wert streben. 

6" bedeutet den ProzeB, die dividierten Differenzen m-ter Ordnung nach den 
n+ 41 Argumenten %», x,,..., x, herzuleiten. Es folgt 











a da  O 0 0 0 0 
ag dag d®aZ 0 0 0 0 
of dat dtat - lat 6" gf 0 0 
lo Th Oh O'W™ORL OH, 6" fo 
fate stots i oe —— Whom = . PP fy aX 
a Joa 8fy08 S%/p08 ... diya gah H+ fal... Hyatt 
Po,1,...,20(%) 1 0 0 ~ © 0 pi 0 (6.0.1) 
a Ja 0 0 0 0 
of dat dta® Mat 0 re 
fo% Sfo% fou ... O" fo % ont, i 0°" fy to 
he — i tas tee O" fo % 3 OT foo tee 5" fy a8 
food Sioa 8/08 ... Moat M+ foak ... Sf, ak 
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Nach einer Laplace-Entwicklung nach den Unterdeterminanten der ersten  Zeilen 
und den algebraischen Komplementen der letzten n+-1 Zeilen des Zahlers, und 
nach den Unterdeterminanten der ersten +1 Zeilen und der algebraischen 
Komplementen der letzten m Zeilen des Nenners, erhalt man 


% da O ... O 00° 2 OF fo ttf, ... PM, 
a2 dag d2af ... O O ... O EFfgaq Mt foay ... P*fgtry 
s=n . ke > : : ° . Ps ‘: 


-o | a oe 
ZI ae bag dah... Saf Stag... Has fy ah? O*+3f, 06-1... Bf, af} 




















at Sat Stat... Slat dtlas ... Mat Sfoak st3fo% ... df al 
K(f, aq, F3fyah,... OT fyah : O41], a8; m) tae 
Im Ausdruck (6.0.2) ist eine neue Schreibweise eingefiihrt worden. Die Deter- 
minante 
UA 2 U3 +++ Am 
U Age og ++. Aom 
W G32 43,3 -++ 43m 
Z Ams Oma +++ Onno 
wird als 
K(u, v, W, ...,2: 4; 5; m) 
bezeichnet. 
Einen knappen Ausdruck leitet man her, wenn man zuerst die Zahler durch 
das Produkt f, /, ... /,, dividiert bevor man die dividierten Differenzen entwickelt. 
Es folgt danach 
Yo,r,..., an (4) = LPR 0 9, 047 Go a8, «=» HF Goad : HFG ab; 2) fot -- fan (6.0.3) 


K(d* +1 fy ag, OT fy ag, ..., HT fy af : Ot fy ax; n) 


mit y,=1//,,7=0,1,...,”. 
Auf eine ahnliche Weise kann gezeigt werden, daB 


n+1 
Po,1,...,2n41(¥%) = > (— 1)*x 


| ag Say usa 0 0 ees 0 
af dag d%a2 ... (e) 0 — 0 
a6 60g taf ... taf dttas ... 8) 
G+} bat? tagtt.. oF Lagt Lt Tagtt ot tantt] K(bfy, fate...» O° feat: IH fait; n +4) 








K (+2 fy ag, H+2 fy a8, ..., MHH2f as: Ht+It1F, af: m) 


(6.0.4) 
und daB 


Po,1,...,2n+1(*) 
— Ko * Poe, MF 70.08, Mt goat: HH G—adin) 4 4 (6.0.5) 
K (+2 f, a, +2 fy ak, ..., M+2fo at: MtIFIA ai: m) (OCR *enth 





gilt. 


— 
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Da nun %o, %4,---) X, 


x,t, "=> (4;)"/n! darf a) durch « = x —¢ ersetzt werden, 
2, (x), das man jetzt mit Fug und Recht als g,,,() bezeichnen kann, 


und p 1...., 
ist durch 
"S (= 1) 4, Bs 
—_ s=0 
Pan(*) = “Sra nl Ce (6.0.6) 





gegeben, wo 














ie 9 0 él 0 0 oe 0 
a2 2a 2! ae 8) 0 nbs 0 
s| a? vor} r(y—1) ar? we W(¥—1) on. (y—s 4+ 2)ar—St1 r(y—1) * (y — s) a"%—s—1 0 
A a” nal n(n —1) a"? ... M(n—1)... (n—s+2)a%—8+1 n(n —1) (n — s) a"—S-2 ind n! wal 
: ITs! 
s=0 
(6.0.7) 
p=n 
— K (pl) pls) — @ #s-2) S\ s—p)(_ 4) (*% — 1 n—1—p. 
K (fff a fi ro DUB (—1) ( ; )p! a 
—s m+j—1) 4n+9-1-9) p | gi-1-?. 6.0.8 
2 | A ts (ore ') pla’ im+1) ) 
rs s) (S)\ #(s—1) (S) 4(s—2) m+J—1) n+j-1). 
K (.() 0. (QA M0: (755 Yani), 
und 
p=n 
= (m+1) (n) m+ 1) sn+1—p) (__ 4)p (% — 1 n—1-p. 
c K (is a—(n+1)f rover Be : )is ( 1)? ( : )plar3-?; 
p=i , 
(nt+i—p) (% +I) (__ 4) (*) 41 2-?- 6.0.9 
Pe : "3 )(—1) (,) Pla nm) ( ) 
=K(t,a,...,08:("FI7 "G41 +4), 
Nach gleichzeitiger Anwendung der Gleichungen (6.0.6), (6.0.7), (6.0.8) und 
(6.0.9) gewinnt man schlieBlich: 
— as fis) ~~ $ fis—)) SH" VS f(S—n) _ 
K > 71 -_ 3 ree Dy a a ic. ') 1)! fri i). n+1) 
(6.0.10) 





an ( )= 
7 K (tas (" FIN Gynt 4) 


4 
Auf eine ahnliche Weise kann man zeigen, daB | 
K(1,0,....0%: ot ‘)@- 1)! ft); m+ 4) 


2n+1 
K (1,04... 0: Paes NG- D1; +4) (fo)?"*7, (6.0.41) 





Fan (x) = 
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und 
S=At] og (sy S=Mt+1 sn s—1) s=nt+l ss xs—m) ; = 
K( SP. y SBR... ERE ret teint 
Ponti(*%) = — _ — 
n 
K(1,0,...,0" C4) eH- +0, :m-+4) 


und schlieBlich, daB 


K (1,0, cong h*s 4 fat ‘) giti-9: a + 2) 


K(t,a,.. ., am: Poa 4). ;m+1) 





Ponti (*) = (fo)°"** (6.0.13) 


4 


ist. Mit den Bezeichnungen 
fs! =4,, (6.0.14) 


Hy = Ay %”, py = a, a"+ aga”, ..., 4, = 4,0" +a,_,0" *+---+a_ (6.0.15) 
sieht man, daB die Gleichungen (6.0.10) und (6.0.12) so geschrieben werden 
kénnen 

Kos Mas +++) Mn? Gi4j—a5 m + 1) 

Pan(2) fon K(a*, a3, ...,13 Gy z_g; @ + 1) 
0 (x) = K (bys Mas +++ nga? Gi4j—-15 2 + 1) 
sats K(a", a®—3, 22, 1: @j45-15" + 1) 





(6.0.16) 





Diese Ausdriicke sind urspriinglich von JAcoBI [20] als die Quotienten zweier 
Polynome n-ter Ordnung (bzw. als die Quotienten eines Polynoms (m-+ 1)-ter 
Ordnung und eines Polynoms n-ter Ordnung) gegeben, deren Potenzreihenent- 


wicklung mit der Potenzreihe >} a,«° bis zum Glied a,,«2" bzw. ag,4,0?"*? 


iibereinstimmt. aes 


Da die Argumente gewissermaBen unendlich nahe beieinanderliegen, wirkt 
der Algorithmus dieser Arbeit wie eine Art Stieltjes-Summation. Durch eine 
ahnliche Rechnung wie bei den Gleichungen (6.0.1) bis (6.0.3) beweist man, wenn 
x,=t+s At, s=0,1,... und At gegen Null strebt, daB gilt 








On). —> (4t)*" 02, (2), (6.0.17) 
oft) —> (At)"**(m +1)! ens ld), (6.0.18) 
mn) (4 t)® ” Moy (t ), (6.0.19) 
ment) _» (At)"**(m +1)! Menai (lt), (6.0.20) 
wn > (At)"n! we, (t), (6.0.21) 
oft —> (At)** 2 a.,41(2), (6.0.22) 

_- K (ay, Ag, «.., Oy? Aj 4 5-15 M) 
Man—1 (0) K(ay, _ Ay: _ m—1)’ (6.0.23) 
die, (f) = K (dq, Gg, ---, @n4y : G47; M) (6.0.24) 


Ki (Gy, Bgs ++ +> Oy Ba g—y3 8) 





(6.0.12) 
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ai ree on a (60.35 
Can er? (6026) 
nll) = eer a a (60.27 
Fans ll) =e Ee (6.0.28) 


Mit Hilfe dieser Ausdriicke kann man die Determinantenausdriicke (4.10.1), 
(4.10.8) (durch Anwendung der Beziehungen (4.1.9) bis (4.2.3)) herleiten. Man 
kann eine groBe Zahl von Differenzen-Differentialbeziehungen zwischen diesen 
GréBen herleiten ; um nur eine zu geben, z.B. 








d {er} — OF (9+ A) yp Wy (6.0.29) 


dt Oy Oy o? 


Aber in gewissem Sinne ist diese Entwicklung der Theorie schon von FROBE- 
nius [21] durchgefiihrt worden. Er definiert 











Coy = K(dg_y 44, Gey as -++1 Og? Mey 4 j45—15%), (6.0.30) 
= K(a,a°+a,_,a°*+--+,a, 0% +a, ga°* +--+)... a, abt ease 
H+ ggg ey 5 pgs +4), 
Sey =K(1,0,...,0"2 deyy-j-4435¥ +1), (6.0.32) 
so daB, z.B. 
Mana (!) =<——, @aq(t) =—Bte, (6.0.33) 
an(t) =z, Says (l) = ——BEE, (6.0.34) 
Tey (t) = — et, TMan+1(t) = — Apia (6.0.35) 
FROBENIUS gab (unter anderen) die folgenden Resultate 
Cert Re_19 — Cey41 Re y1 = Ce, Re, (6.0.36) 
Cevtt Resa» — Cosa Re 41 = legit eRe, - (6.0.37) 
Cetiy Rey — Ce» Rerry = Cetirti eRe ya, (6.0.38) 
Cevta Rey — Coy Rey4t = Cetry410¢Re_1 5, (6.0.39) 
Risege (Se = —(F+9+1) Cop 4 Re vga, (6.0.40) 
Ce, Lisear we —(E+%+1) 644542 Re-1,, (6.0.41) 
chase gp—<Btp er = (E+ 9+ 1) eee, (6.0.42) 


Cé+1,¥% 
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c2 d Re, ye (€ ++ 1) Cey Re y41 
&,»t+l Gy Coop garrett of t*t+2 ’ 








(6.0.43) 


d d 
Ceti» gy evta — %Ev+a Ge +1,» = (E+¢+1) Ce Cera v4, (6.0.44) 


d Ce+ ” d Ce, ae 
Cy a TTT ~ 041 GE Epp ~ -e EH ety (6.0.45) 


Cc d e+ Cc @ ie == ¢ Cc 
&,y at Etyvt+1 &,v+1 dt Et» €—-1,9%§+1,941° 





(6.0.46) 


Durch Anwendung der Gleichungen (6.0.36) bis (6.0.46) kann man ein Differential- 
Analogon zu den Gleichungen (4.1.9) bis (4.2.3) herleiten. Zunachst bildet man 
die Folge der Determinanten c¢,+, ,, ¢,41,,41, 7 =0, 1,... durch Anwendung der 
Gleichungen 


“4 d C,, ad Cy 
Cr+iyr (C,-1,¢—-1) . {67-1 at — va ary sa 7 = 0, 1, eeey (6.0.47) 


ad ad Criiy d Cry 
Cr+ijrt+i— (c, »—-1) : {e,,+ adi err — Cot ay dt 3s} r=0,1,..., (6.0.48) 


und danach die Polynomfolgen S,,, S,,41, 7,,,, 7,41, durch Anwendung der 


Gleichungen 
7+1 2 
R os a2rt+ Cry d Ry 9—1 = O, a, core (6.0.49) 


100 BP lp, gag Ot by ” 





a2rt+2 Gate d R,,+ 


R, +41 - (27+ 1) Cry dt Cr+i,r a27t+l r=0,1,... (6.0.50) 





mit R,,=S,, und 7, ,. SchlieBlich bildet man die GréBen Pon(%), Ponsa (%) in 
den Gleichungen (6.0.16) durch Division. 


Es sollte jedoch erwahnt werden, daB es viel weniger umstiandlich ist, zunachst 
die aufeinanderfolgenden Ableitungen /) (t)/s!, s =0, 1, ... herzuleiten, un. 7anach 
durch Anwendung des g — d-Algorithmus oder einer anderen geeigneten Methode 
die Kettenbruchentwicklung von /(x) zu gewinnen. Die rationalen Quotienten 
(6.0.16) sind dann die u-ten bzw. die (n+ 1)-ten Naherungsbriiche dieses Ketten- 
bruchs. 
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Uber die numerische Behandlung von 
Integralgleichungen nach der Quadraturformelmethode 


Von 


HELMUT BRAKHAGE 


Einleitung * 


Die Arbeit befaBt sich mit der numerischen Behandlung von nichtsingularen 
Integralgleichungen zweiter Art nach der Quadraturformelmethode. Nach einigen 
einfiihrenden Betrachtungen in Nr. 1 wird in Nr. 2 eine neue Art der iterativen 
Korrekturrechnung angegeben, die fiir numerische Zwecke besonders geeignet 
erscheint. Davon unabhangig werden in Nr. 3 Fehlerabschatzungen hergeleitet. 
Die letzte Nummer enthalt Bemerkungen zur praktischen Durchfiithrung mit 
Beispielen. 

Die Korrekturrechnung ergibt sich im Zusammenhang mit dem Ubergang 
zu einer verfeinerten Quadraturformel mit erhéhter, etwa verdoppelter Stiitz- 
stellenzahl. Verglichen mit dem iiblichen Stérungsansatz, wie er von H. BUCKNER 
fiir die Methode des Ersatzkernes angegeben wurde ([J], S. 95), wird hier vor 
dem eigentlichen Korrekturschritt jeweils noch ein vorbereitender Schritt er- 
forderlich, der eine Glattung der Defektfunktion bewirkt. P. A. SAMUELSON hat 
in einem allgemeineren Rahmen die iterative Lésung von Integralgleichungen 
unter Benutzung einer Naherung fiir die Resolvente behandelt [7] und dabei 
auch die Anwendbarkeit auf Stiitzstellenmethoden angedeutet. Hier wird die 
explizite Aufstellung einer Resolventennaherung vermieden. Die Fehlerabschat- 
zungen verallgemeinern einen Ansatz von L. W. KANTOROWITSCH [4]. Beziiglich 
weiterer Literaturhinweise, insbesondere zur Fehlerabschatzung, sei auf [1], [4] 
und [5] verwiesen. Was die Darstellung anbetrifft, so wird fiir die nach der 
Quadraturformelmethode gewonnenen Niaherungsgleichungen eine Operator- 
schreibweise eingefiihrt, die fiir die theoretische Behandlung vorteilhaft ist. 


Auch auf die Bestimmung von Eigenwerten und Eigenfunktionen von Inte- 
gralgleichungen nach der Quadraturformelmethode lassen sich die benutzten 
Uberlegungen iibertragen, worauf in einer gesonderten Arbeit eingegangen 
werden soll. 

1. Problemstellung 


Zugrunde liege die Integralgleichung 


(4) x (s) +f R(s,t) x(t) dt — f(s) =0 (0SsS1). 
0 





* Fiir verschiedene Anregungen und kritische Bemerkungen bin ich Herrn Prof. 
Dr. J. WEISSINGER und Herrn Priv.-Doz. Dr. W. WALTER zu Dank verpflichtet. 
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Indem man das Integral von (1) naéherungsweise durch eine Quadraturformel 


(2) fuQdtw 3 A,u(t) 


0 j=1 


auswertet, wird (1) in ein lineares Gleichungssystem 
(3) y (4) + LAH) A; y(t) —f(4)=0 (¢=1,...,m) 
17= 


iibergefiihrt. Fiir das Folgende ist es zweckmaBig, mit Hilfe der Stiitzstellen- 
werte y(¢;) nach einem Vorschlag von E. J. NystR6m [6] durch die Festsetzung 


(4) y(s) =f(s) — x k(s, 4) Ay y(t) 


eine fiir alle Argumente s€ [0, 1] definierte Naherungslésung einzufiihren. 


Will man von vorneherein sichern, daB die Naherungslésung (s) die exakte 
Lésung x(s) mit hinreichender Genauigkeit approximiert, so wird man haufig 
die Stiitzstellenzahl m unnétig hoch wahlen miissen. Da der Rechenaufwand 
fiir die Auflésung von (3) mit steigender Stiitzstellenzahl sehr rasch ansteigt, 
wird es oft ratsam sein, zundchst versuchsweise ein Gleichungssystem mit gerin- 
gerer (etwa halber) Stiitzstellenzahl m’ aufzulésen, die Genauigkeit dann aber 
dadurch zu iiberpriifen, daB man eine Einsetzprobe in das Gleichungssystem (3) 
mit der genaueren Quadraturformel macht. Ist das Resultat der Einsetzprobe 
nicht befriedigend, so kann man — anstatt dann das Gleichungssystem mit 
der héheren Stiitzstellenzahl direkt aufzulésen — die erste Naherung unter 
Benutzung der bei der Einsetzprobe sich ergebenden Defekte durch eine einfache 
Korrekturrechnung verbessern. Es soll im Folgenden angegeben werden, wie eine 
solche Korrekturrechnung zweckmaBig durchgefiihrt werden kann. Vor dem 
eigentlichen Korrekturschritt erweist sich ein vorbereitender Schritt als erforder- 
lich, der eine Defektglattung bewirkt. Bei der Korrekturrechnung werden dann 
in das Gleichungssystem mit der geringeren Stiitzstellenzahl lediglich neue in- 
homogene Glieder eingesetzt. Es ergibt sich kein wesentlicher neuer Arbeits- 
aufwand, da die Koeffizientenmatrix von der friiheren Auflésung her etwa als 
Produkt zweier Dreiecksmatrizen vorliegt. Die Korrekturrechnung kann als 
erster Schritt eines i.a. sehr rasch gegen y(s) konvergierenden Iterationsver- 
fahrens gedeutet werden, dessen Konvergenzeigenschaften abgeschatzt werden 
kénnen. Oft wird der erste Korrekturschritt ausreichen. 


Es sei 


1 m’ 
(5) Ju dtw 2 Aj ul) 


j=1 


die grébere Quadraturformel mit der geringeren Stiitzstellenzahl m’, und 


(6) y(t) + DAG.) ASG) — HG) = 0 = Ay.) 
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das zugehérige Gleichungssystem. Entsprechend zu (4) sei y’(s) definiert durch 
a , ‘ , 
(7) y'(s) = f(s) — Lk (s, 4) Aj ¥'(G)- 


Die Darstellung des Folgenden wird in Operatorschreibweise iibersichtlicher. Durch 
den Kern k(s,¢) wird ein Integraloperator K in der tiblichen Weise definiert: 


(8) Ku=f k(s,t) u(t) dt. 


0 
Dann ist (1) gleichwertig mit 


(9) (I+ K)x—f=0 


(I identischer Operator). Die den Operator K approximierenden Niaherungs- 
operatoren H, H’ seien definiert durch 


m m’ 
(10) Hu=)DR(s,t;) Aju(t), H'u=> k(s,t) Aj u(t). 

j=1 j=1 
Die naherungsweise Auflésung von (9) nach (3) und (4) bzw. (6) und (7) ist 
offenbar gleichwertig mit der Auflésung der Operatorgleichungen 


(11) (i+ H)y—f=0 
bzw. 
(12) (I+ H’) y’—f=0. 


Die Gleichung (11) ist namlich identisch mit (4), und diese wiederum enthalt 
(3) fir s=¢; (t=1,..., m). 

Der Einfachheit halber werde vorausgesetzt, daB k(s,¢), f(s) und alle auf- 
tretenden Funktionen stetig seien; doch gelten die meisten Ausfiithrungen auch 
allgemeiner, sofern nur die Anwendung der Quadraturformelmethode noch sinn- 
voll ist. Weiter wird angenommen, daB die Gleichungen (9), (11), (12) eindeutig 
auflésbar sind, d.h. daB die Umkehroperatoren (J+ K)+, (J+ H)7, ([+H’)7 
existieren. SchlieBlich werden noch die folgenden Normen 


|||] = Max {|«(¢)||oS¢< 1} 
|7 |] = sup {|Z el | Il] = 1} 
(T linearer Operator) eingefiihrt. 


(13) 


2. Iterative Bestimmung einer Naherung 


Es sei y die zu approximierende Lésung von (11), yo eine beliebig vorge- 
gebene Naherungslésung, die als Ausgangsnaherung fiir die Iteration dienen soll. 
Fiir yp kann etwa die Lésung y’ der Gleichung (12) gewahlt werden. Der sich 
beim Einsetzen von yp in die Gleichung (11) ergebende Defekt werde mit 7, 
bezeichnet, 


(14) 1 = (1+ H) ¥—f. 


Numer, Math. Bd. 2 14 
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Setzt man fiir f nach (11) (J+) y ein, so folgt fiir die zu bestimmende Korrektur 
ye 
r= (1+H)(y—y), dh. y—¥=—(I+H)*1. 


Nach dem iiblichen Vorgehen der Korrekturrechnung ware es naheliegend, y — v9 
zu approximieren durch — (I+ H’)+7,. Die Werte dieser Naherung fiir y — yy 
an den Stiitzstellen ¢; geniigen einem Gleichungssystem vom Typ (6), in das als 
inhomogene Glieder nur die Stiitzstellenwerte 7,(t;) eingehen. Man kann hieraus 
schlieBen, daB — (I+ H’)+7, nur dann eine brauchbare Naherung fiir y— yo 
ergeben wird, wenn die Defektfunktion 79(s) einen hinreichend glatten Verlauf 
hat. In der Praxis zeigt sich, daB die Defektfunktionen im allgemeinen stark 
oszillieren, so daB eine Modifikation des angegebenen Ansatzes fiir die Korrektur- 
rechnung erforderlich wird. 


In einem vorbereitendem Schritt wird y) durch 


vo =Vo—%o 


ersetzt. Die Naherung yg wird im allgemeinen nicht besser sein als yy. Be- 
trachten wir jedoch den Defekt 7 der Funktion yf, so gilt 


rm = (1+ A) ¥¢ —f=(1+ A)» —f—(+ A) mn =—H; 


d.h. durch den Ubergang zu yj wird eine (durch den Operator H bewirkte) 
Glattung des Defektes erreicht *. 


Ersetzt man jetzt in der Gleichung 

y¥—¥ =—(1+H)* 

den Operator H durch H’, so erhalt man die verbesserte Naherungslésung 
n= — (+H) 

oder, ausfiihrlich geschrieben 

(15) w=Y—%+(Ul+H')*Hr mit m=(1+H)y—f. 

Die iterative Anwendung dieses Vorgehens fiihrt auf die Iterationsvorschrift 

(16) Ynt1=In—Mm+(I+H')* Hr, mit 7,=(I+H) y,—f. 


Die nach diesem Verfahren berechneten Naherungen y, konvergieren im all- 
gemeinen sehr rasch gegen y, wie im Folgenden gezeigt werden soll. 


Die durch (16) definierte Funktionenfolge y,, konvergiert, falls 
(17) || + 4’) (H’ — A) H|| <1 
ist, gegen die Lésung y der Gleichung (11). 
Ist a= > 14jl und gelten die Abschatzungen 
jen 


(18) |(K —H),k(r,t)| <6, |(K —H’),k(t,t)| <8, 
* Vgl. hierzu die beiden Beispiele in Nr. 4. 
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so gilt 
(19) |(H’ — H) H||<a«(6+6’). 


Wenn die Quadraturformel (2) fiir «const. exakt ist und die A;=0 sind 
— diese Voraussetzungen sind bei den gebrauchlichen Quadraturformeln er- 
fiillt — so ist xa=1. Die in (18) auftretenden GrédBen sind nichts anderes als 
Quadraturfehler; es ist z.B. 


(K — H), k(t, t) wells t) k(t,t)dt — 3A; k(s,4) (tt). 
F - 


j=1 


Der Forderung (17) kann demnach durch hinreichend genaue Quadraturformeln 
Geniige getan werden. 


Zum Beweis der Behauptungen multiplizieren wir die Gleichung (16) mit 
I+H’. Nach einfacher Umformung ergibt sich 


(20) (I+ 1’) \ngs =f + (H'— A) f—(H’— A) yy. 
Unter der Bedingung (17) hat also die Gleichung 
(21) (I+ H’)z=f+(A’— HA)f—(A’—#A)Hz 


genau eine Lésung z, gegen die die Iterationsfolge y,, konvergiert. Nun ist abér 
(1+ H) y=f, und deshalb 


_ (+) y=f4+(H'— A) y. 


Fiihrt man auf der rechten Seite die Substitution y—/— Hy durch, so sieht 
man, daB auch y die Gleichung (21) erfiillt, d.h. daB z=y ist. Die Abschatzung 
(19) leitet sich unmittelbar ab aus der Formel 


(H’ — H) Hu=(H' — H) 


ime 


k(t, t;) A; u(%) 


t 
1 


= [(K — H) —(K —H’)]-R(t, 4) Aju(t). 


j=1 


Es sei noch angemerkt, daB unser Iterationsverfahren (16) auch als verein- 
fachtes Newtonsches Verfahren zur Lésung der Gleichung (11) interpretiert 
werden kann. Die vereinfachte Newtonsche Iterationsvorschrift lautet 


Vn+1 _ Vn = Gr,,, 


wobei G ein Naherungsoperator fiir (+H) ist. Man erhalt gerade die Vor- 
schrift (16), wenn man 


(I+ H)t=I-—(1+H)1H 
durch 
G=I-—-(I1+A’')*H 


Numer. Math, Bd. 2 14a 
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approximiert. Noch kirzer gesagt: Die Ersetzung von H durch H’ wird nicht 
in (I+ H)+, sondern bei der Resolvente vorgenommen. 

Auf die praktische Durchfiihrung der Iteration wird in Nr. 4 eingegangen. 
Sie ist einfach und fiir die Programmierung auf elektronischen Rechenmaschinen 
geeignet. In Nr. 4 finden sich auch Angaben tiber die Aufwandsersparnis bei 
dem beschriebenen Vorgehen. 


3. Fehlerabschatzung 


Die hier durchgefiihrten Fehlerabschatzungen sind unabhangig von den Aus- 
fiihrungen der vorangehenden Nummer. Es sei eine Naherungslésung %, fiir die 
Integralgleichung (9) vorgegeben. Wir nehmen wieder an, daB (9) genau eine 
Lésung x besitzt, d.h. daB (I+ KK) existiert. Theoretisch ist es einfach, mit 
Hilfe des Defektes beziiglich der Gleichung (9) 


(22) Rx, = (1+ K)x,—f 

die Differenz x — x, abzuschatzen; es gilt namlich 

(23) (I + K) (% — %) = —R%, 
also 

(24) llx — %oll S||(Z + K)* |] - [1B xl]. 


Damit ist die Frage der Fehlerabschatzung darauf zuriickgefiihrt, eine Schranke 


(25) (7+ K)4|| Sop, 


sowie eine Abschatzung fiir Rx, zu bestimmen. Wir beschaftigen uns zundchst 
mit der zweiten Frage. Die Berechnung von Rx, erfordert eine im allgemeinen 
nicht explizit ausfiihrbare Integration. Deshalb wird eine Quadraturformel be- 
nutzt — sie sei durch (2) gegeben — und der Defekt 7, beziiglich der entsprechen- 
den approximierenden Gleichung (11) berechnet: 


(26) 19 = (I+ H) %—f. 
Es gilt dann die Formel 
(27) R %)=%+(K — A) %, 


die jedoch zur Abschatzung noch nicht gut geeignet ist. Man kann durch eine 
Modifikation dieser Betrachtungen zu einer wesentlich besseren Abschatzungs- 
formel gelangen. Betrachtet man namlich statt x» als HilfsgréBe die Naherung 
(28) x = %—% 

(die Einfiihrung von x¢ ist, nebenbei bemerkt, véllig analog zu der von y¥ in 
der vorangehenden Nummer), so gilt entsprechend zu (26) und (27) 

(29) R xg =7$ + (K — A) xf, 

wobei 

(30) rm = (I + H) xf —f 
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ist. Wegen ; 
rf =(I +H) %—f—(I+H)m=—H1 
und 
xe = %—%=/—H X% 


ist (29) Aquivalent mit 
(34) Rx =—Hr+(K—A)f—(K—HA)A x. 
Diese Gleichung bildet die Grundlage fiir die folgenden Abschatzungen. 
Ist wie in Nr. 2 

(32) a=SIAi], (KH) A(e 01S, 
so folgt aus 

\|(K — H)Hu| = 3 (K — H), k(t, t;) A; u(é;) 
die Ungleichung = 

|(K—H)Hu|s a 6 Max'|(i)|. 


Hat man dariiber hinaus die Ungleichungen 
(33) (K—H)f|<e, Max{ }|4;k(s,4)||o<ssib <M, 
j=1 


so ergibt sich aus (31) die im ganzen Intervall [0, 1] giiltige Abschatzung 
(34) |R xg | SM Max|79(4)| + € + a 6 Max |x9(t; )|. 


Wenden wir nun die (fiir jede Naherung giiltige) Formel (24) auf xf an, so folgt 
schlieBlich die gewiinschte Abschatzungsformel 


(35) | x(t) — x(t) + 79(t) )| SP [M Max |79(t) | + & +a 6 Max |x9(4) (t;)|] 
mit 
7o(t) = (I + A) x — f= xo (t) )+ Day ;) A; X(t) — f(t). 
7= 
Bringt man in (35) 79(¢) auf die rechte Seite und beschrankt sich auf die Ab- 
schatzung an den Stiitzstellen, so erhalt man die einfachere Formel 


(36) |x(t;) — x(4)| SP (“+ ; ) Max |70(¢; )| + +a 6 Max | Xo (t,) |. 


Ist x die Lésung und x, eine Naherungslisung der Integralgleichung (9), so 
gelten die Fehlerabschatzungen (35) und (36). Eine Schranke p fiir die Norm 
des inversen Operators (I+ K)+ ist im Falle ||K||<1 durch p=(1—||K||)+ ge- 
geben. Ist diese Bedingung nicht erfillt und ist die zum Gleichungssystem (3) 
gehérige inverse Matrix 


B= {bj} = {8,4 — (tir) A}? 
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bekannt, so kann man (falls der folgende Nenner positiv ist) 





— 1+ ||AI]4+6 M) ra * 
(37) p= 1— (1+8 M) 5 mut B = Max 2/4, 
setzen. Die GroBen a, 6, &, M sind in (32), (33) erklart. 


Es sei darauf hingewiesen, daB diese Abschatzung der Norm von (J+ K)* 
ganz unabhangig von den Formeln (35), (36) gilt. Das bedeutet, daB die in (35), 
(36) den GréBen 7,, «, 6, e, M zugrunde liegende Quadraturformel nicht identisch 
zu sein braucht mit der Quadraturformel, auf die sich B, B, 6, M in (37) beziehen. 
Man kann sich z.B., wenn man eine Naherung nach dem Verfahren von Nr. 2 
bestimmt, bei der Fehlerabschatzung (35) auf die Quadraturformel (2), bei der 
Abschatzung (37) auf die Quadraturformel (5) stiitzen (letzteres ist vorteilhaft, 
weil in diesem Falle die der Quadraturformel (5) entsprechende Matrix B als 
bekannt angesehen werden kann). 


Der Beweis fiir (37) geht aus von der Identitat 
(I +H) —K](+ K)=([+ 4) —(K—A#A)K, 


aus welcher 

(38) (i+ K)*= [I — (1+ 4)*(K — A) K}7 (I — (I+ H)*K) 
folgt. Nun ist, wie wir sogleich zeigen werden, 

(39) (7+ H)3||<1+pM. 


Ist namlich w= (J+) v, so haben wir einerseits 


v (4) = Sdn mlt), also |v(t)|<B||x||, 


andererseits 


also 
lle | S I] || + M Max |o(¢)| S (1+ BM) ||4|]. 


Aus der Gleichung (39), die iibrigens in 4hnlicher Weise bei KANTOROWITSCH [4], 
S. 101 benutzt wird, folgt unmittelbar die Ungleichung 


|Z + H) 7 (K — H) K||/S(1+6M)6 


und hieraus mit Hilfe von (38) die Gleichung (37). 


SchlieBlich sei noch ein Weg skizziert, wie man zu Schranken e, 6 gelangen 
kann. Nimmt man an, da8 der Quadraturfehler 


1 m 
E(u) a u(t) dt — 2,4; (4) 
j= 
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fiir eine /-mal stetig differenzierbare Funktion u die Abschatzung 
(40) |E(u)| SC Max | u (z)| 


gestattet, so folgt aus der Regel von LEIBNIZ 
l 
(41) | E(u, + u)|<C Dt ) Max | u{? (z)| Max|uf-?)(zx)|. 


Die beiden GréBen (K —H)f, (K —H),k(t, t) sind nichts anderes als die Qua- 
draturfehler der Funktionen 


u(t) =k(s,t) f(t) bzw. u(t) =R(s,7) R(t, 2), 


die mit Hilfe von (41) abgeschatzt werden kénnen, wenn / und k /-mal stetig 
differenzierbar sind. 


Wenn man # nach (37) und die GréBen ¢ und 6 in der angegebenen Weise 
bestimmt, so ist unsere Abschatzungsformel (35) im Sonderfall 7,=0 mit der 
von KANTOROWITSCH [4], S. 102 gegebenen etwa gleichwertig. Sie hat aber den 
Vorteil, auch fiir solche Naherungen x) anwendbar zu sein, welche nicht durch 
exakte Auflésung eines linearen Gleichungssystems gefunden werden. So wird 
man, wenn ||H||<1 ist, z.B. durch Berechnung einiger Glieder der Neumann- 
schen Reihe fiir (1-++ H) + zu Naherungslésungen gelangen, welche im allgemeinen 
nicht die Bedingung 7,=0 erfiillen. Oder man wird mit Hilfe des in Nr. 2 ge- 
schilderten Verfahrens Naherungslésungen y, erhalten, fiir die 7, (das entspricht 
7%) nicht verschwindet. 


Weiter bietet die von der speziellen Fehlerabschatzung (40) unabhangige 
Definition von 6 und « in vielen Fallen Vorteile. In diesem Zusammenhang sei 
auf die von Davis [2], [3] fiir analytische Integranden angegebenen Abschiat- 
zungen des Quadraturfehlers hingewiesen, welche keine Ableitungen enthalten 
und deshalb fiir unseren Fall, wo der Integrand ein Produkt ist, besonders ge- 
eignet sind (vgl. dazu das zweite Beispiel zur Fehlerabschatzung in Nr. 4). 

Es sei noch darauf hingewiesen, daB zur Abschatzung von x— x» an den 
Stiitzstellen gemaB (36) die Kenntnis von 7, lediglich an den Stiitzstellen vor- 
auszusetzen ist. 


4. Bemerkungen zur Korrekturrechnung und Beispiele 


Bei dem Verfahren von Nr. 2 wird man die Stiitzstellensysteme ¢;, ¢; nach 
Moglichkeit so wahlen, daB die ¢; unter den ¢; enthalten sind; doch ist dies z.B. 
bei der Verwendung GauBscher Quadraturformeln nicht erreichbar. 

Nachfolgend wird eine ausfiihrliche Beschreibung der einzelnen Schritte fiir 
die numerische Durchfiihrung der Korrekturrechnung nach (16) gegeben. Es 
sei y,(s) fir s=é¢, (i=1,...,m) und s=tj (i=1,..., m’) bekannt. Dann ergibt 
sich 7,,(s) durch Einsetzprobe mit der verfeinerten Quadraturformel. 


(42) 47,(s)= 8) + Dass )A;¥,(¢;)) —f(s) fiir s=t; und s=¢;. 
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Fallt dabei 7,(s) so groB aus, daB ein neuer Iterationsschritt gerechnet werden 
muB, so besteht die eigentliche Korrekturrechuung darin, die Funktion 


z,(s) =([+ H')7 Hr, 


zu bestimmen, d.h. das Gleichungssystem 


Mm 
(43) 2» (tj) + DR tj) Aj 2n (Gj) = (G4) (¢=1,..., m’) 
I= 
mit der rechten Seite 
(44) ra(s) = Hr, = DR(s,t;) A;7,(t;) fiir s=t; 
j=1 


aufzulésen. Daraus ergibt sich dann y,.4,(s) gemaB der mit (16) gleichwertigen 
Formel 


(45) Vn+1 (s) ™ Fy (s) “te (s) + Zn (s) 


an den Stellen s=¢;. Benétigt man auch die Werte y,,,,(¢;) (etwa um eine neue 
Einsetzprobe zu machen und das Iterationsverfahren fortzusetzen), so hat man 
noch 7#(s) fiir s=¢; nach (44) und sodann 


(46) z,(s) =7%(s) — > k(s, t;) A; z,(%) fir s=¢, 
jul 


zu berechnen. Diese Bestimmung der z, (¢;) erfordert einen verhaltnismaBig groBen 
Rechenaufwand. Es ist wesentlich einfacher und fiir praktische Zwecke meist 
ausreichend, die z,, (¢;) statt dessen durch Interpolation aus den z, (¢;) zu gewinnen*. 
Die Funktion z,(s) hat némlich im allgemeinen ebenso wie r*(s) einen glatten 
Verlauf; der Fehler durch die Interpolation wird daher klein gegen z,, sein. 


Im Normalfall wird man sich mit einem Korrekturschritt zufrieden geben 
kénnen. Die ganze Rechnung setzt sich dann aus folgenden Schritten zusammen: 
Zuerst wird die Gleichung (6) aufgelést und die Lésung y’ als Ausgangsnaherung 
Yo gewahlt. Dann wird y’=y, nach (7) an den Stellen ¢; berechnet und danach 
die Einsetzprobe gemaB (42) vorgenommen. Ist das Ergebnis der Einsetzprobe 
unbefriedigend, so werden die GréBen 7% (t,), z9(¢;) nach (43), (44) und hieraus 
y, (t;) nach (45) bestimmt. 


Der Rechenaufwand fiir die direkte Auflé6sung von (3) soll mit dem Aufwand 
fiir den eben beschriebenen Normalfall unter der Annahme m’ = m/2 iiberschlagig 
verglichen werden. Es wird vorausgesetzt, daB die ¢; unter den ¢; enthalten sind. 
Als MaB fiir den Rechenaufwand wird die Anzahl A der Multiplikationen gewahlt, 
fiir die Auflésung eines Gleichungssystems von » Unbekannten wird A =n°3/3 
gesetzt. Fiir die direkte Auflésung von (3) mit der feineren Unterteilung hat 
man A,=m*/3, fiir die Bestimmung der yg (¢;) nach (6) mit der gréberen Unter- 
teilung A,—m'/24. Die Einsetzprobe (42) (einschlieBlich der Berechnung der 





* Dagegen ware es im allgemeinen nicht sinnvoll, die Interpolation direkt bei 
den Werten y,,(¢;) vorzunehmen. 
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Yo(t,)) erfordert A,=5m?/4, die Korrekturrechnung A,=3m?/4 Multiplikationen 
(man beachte, daB die Auflésung von (43) nur A = m?/4 erfordert, da die Koeffi- 
zienten von (6) und (43) identisch sind). Bestimmt man noch die 29(¢,) und 
y, (¢;) in der strengen Form nach (46), so hat man fiir A, den gréBeren Aufwand 
A,=7m?/4; benutzt man fiir die z)(¢;) quadratische Interpolation, so erhéht 
sich die Anzahl der Multiplikationen nur unwesentlich um 3m/2. Es ergibt sich 
mit dem erstgenannten Wert von A, 

1 


mM ay — 


| 


A,:4,=— mw'im: A,:4,= 


_ yy — mM. 
15 4+ 


\o|+ 


Im folgenden 2. Beispiel (m=60) ist also A,;:A,~%15:1, A,: A, 30:1. 
1. Beispiel. Die Integralgleichung 
1 
x(s) + s f _*(di_ _y, 


a 1+(s—t)? 
-1 


ergibt unter Beachtung von x(— ¢t) = x(t) 





1 
(47) “(+5 { Eeaom + aera AO dt=t. 
0 

Die Anwendung der Simpsonregel mit den Stiitzstellen ¢;: 0, $, 1 und die Auf- 
lésung der entsprechenden Gleichungen (6) fiihrt zu den drei Naéherungswerten 
Yo (t;) (s. unten). Als verfeinerte Quadraturformel verwenden wir wieder die 
Simpsonregel mit ¢;: 0, %, +, =, 1. Die Ergebnisse der Einsetzprobe und der 
Korrekturrechnung sind in der folgenden Tabelle zusammengefaBt. Die GréBe 
des Defektes 7, zeigt, daB y, praktisch eine Lésung von (3) ist. 


























t Yo 10*+ 10* + rF 10*+ z5 V1 10’-7, 
0 0,658 3092 9,130 0,286 0,210 0,6574172 4 
4 0,663 9298 1,202 0,260 0,186 0,663 8282 4 
$ 0,682438 3 — 7,186 0,209 0,140 0,683 1709 2 
q 0,7148314 —0,258 0,177 0,116 0,7148688 2 
1 0,7565016 1448 0,166 0,114 0,7577358 1 


Fiir die Fehlerabschaétzung bemerken wir, daB 








aie 4 : _ ; 
“rT Im(+) mit z= t+1n, 
also 
dxvA 4 = A! ' 
“dh 44%? = |ImGar),-1 Si! 
und deshalb 
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ist. Ferner ist 
1 1 
\|K || = Max f|k(s,#)| dt = f k(0,t) dt=2, 
s 0 0 


also p= (1— ||K||)?=2. 

Wir fiihren eine Fehlerabschatzung in drei Stufen durch. Sind nur die drei 
Naherungswerte y (¢;) bekannt, so gilt fiir die verwendete Simpsonregel die 
Abschatzung (40) mit /=4, C= 1/2880. Bezogen auf aiese grébere Quadratur- 
formel ergibt sich gemaB (41) mit (48) 


Sa —- & 7 = or 4 
2880 2 => 0,01 , E 88 0,005 





und daraus nach (35) wegen |¥o(¢;)| 0,76 und 7y,=0 
(49) |x (t) — yo (t)| < 2(0,0054 + 0,017 - 0,76) < 0,037. 


Nun nehmen wir an, die Einsetzprobe sei durchgefihrt, d.h. die beiden ersten 
Spalten der obigen Tabelle seien bekannt. Beziiglich der Simpsonregel mit 
halbierter Schrittweite ergeben sich neue Schranken ¢*=e¢/16, 6*=6/16. Mit 
| 79 (4;)| S 0,001, |y¥o(¢,)|S0,76 und M=2/x folgt nach (36) die genauere Abschat- 
zung 


(50) |x (t;) — Yo (4)| S 0,005. 


Dieses Beispiel zeigt, daB man allein mit der Einsetzprobe die Fehlerabschatzung 
unter Umstanden wesentlich verbessern kann. Hat man schlieBlich alle Werte 
der Tabelle errechnet, so ergibt sich aus (36) entsprechend 


(54) |x (t;) — y, (t,)| S 0,0024. 


2. Beispiel. Fir die Ellipse soll das Dirichletsche Randwertproblem gelést 
werden. Setzt man die Lésung als Potential einer Dipolbelegung der Randkurve 
an, so ergibt sich die Belegungsdichte x(t) aus der Integralgleichung 


1 
-— 4 _ a—b 
(52) x(s) + f papeesacry dt=Sel), y= ty 
0 


(a, b Halbachsen). Es wird der Fall a=9, b=1 (y=4/5) und g(t)=2 behandelt. 
Die exakte Lésung ist x (#)=1/2. Als grébere Quadraturformel wird die Rechteck- 
regel mit m’=30 Stiitzstellen gewahlt (Aj=1/30, #;=7/30, i=1,..., 30). Es 
ergeben sich die Naherungswerte 








Vo (t;) = 0,4993810 fiir i=41,..., 30. 


Geht man — wieder mit der Rechteckregel — zur doppelten Stiitzstellenzahl 
iiber (A;= 1/60, ¢;=1/60, i=1, ..., 60), so ergeben Einsetzprobe und Korrektur- 
rechnung die Resultate 

Yo(t;) = 0,5006175 — 0,001 2365 (— 1) 
ro (t;) = 0,001 2365 — 0,001 2395 (— 1) 
rp (t;) = 0,001 2365 — 0,000003 0(— 1)* 
(— 4) 
(— 4) 


i 


t 


29 (¢;) = 0,0006205 — 0,000004 5 i 
4, (t;) = 0,500001 5 — 0,000001 5 


1 
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Die einzelnen Schritte der obigen Rechnung lassen sich unter Benutzung der 
Kerndarstellung 


k(s,t) =1+2 > y’cos22v(s + 2) 
v=1 


explizit ausfiihren. 

Dieser Fall einer sehr langgestreckten Ellipse mit dem Achsenverhiltnis 9:1 
ist fiir die numerische Behandlung ausgesprochen ungiinstig. Die obigen Resultate 
zeigen, daB der Fehler durch die Korrekturrechnung etwa um den Faktor 500 
kleiner geworden ist. Eine Fehlerabschatzung ergibt hier nur sehr grobe Schran- 
ken; sie wird deshalb nicht wiedergegeben. 

Zur Durchfiihrung einer Fehlerabschatzung kniipfen wir an die Integral- 
gleichung (52) an und wahlen darin 

a=2, b=1, y=%, 8(s)=76(5+3cos4zs). 

Dieses Beispiel wurde in der Literatur mehrfach nach der Quadraturformel- 
methode behandelt (NysTROM [6], KANTOROWITSCH [4]). Eine Fehlerabschatzung 
fiir die Quadraturformelmethode wurde bisher nicht durchgefiihrt. Es soll nun 
ausgefiihrt werden, daB die Ungleichung (35) in Verbindung mit Abschatzungen 
fiir den Quadraturfehler von Davis [2], [3] zu brauchbaren Fehlerschranken 
fiihrt, wobei die miihsame Abschatzung der Ableitungen von k und g vermieden 
wird. 

Ist u(z), z=t-+7y, reell fiir reelle z, periodisch mit der Periode 1 und ana- 
lytisch in dem Streifen |7|<A, so gilt fiir den Quadraturfehler 

1 
ons _ A m j 
E(u) = fue t)dr— + i (2 Z 
; = 
bei der Rechteckregel mit m Stiitzstellen nach Davis [3] 


(53) | E(u)| S 2e-274™ (4 — e328 Am) 1 Max |u(z+7A)|. 


Nun ist wegen | cos 27(t+7A)|<cosh 227A 





r 1—y* 1+? 
<= 4 it 
|k(s,t+2A)|S ers YTS falls cosh2aAS - 


ist, und 


le(t+2A4)| See (5 +3 cosh2zA). 


Legen wir unserer Fehlerabschatzung die auch in [4] gewahlte Stiitzstellenzahl 
m= 16 zugrunde, so ergibt sich 


6=2,8-10°, e=2,1-1077; 


dabei wurden fiir 6 und ¢ in (53) verschiedene Werte A benutzt, die sich aus 
afi 2 _— 
e st (1+ Z)y bzw. e~?74— =(1+2 m)Y 


ergeben. Unter der Annahme, daB x(t) die exakte Lésung des entsprechenden 
linearen Gleichungssystems ist, hat man dann die Fehlerschranke 

(54) |x(t) — xo(t)| S$ p(2,1- 10-7 + 2,8-10°5-0,1) S p- 3,1 -10°S 

(es ist namlich |x (¢,)| 0,1). 


196 Berichtigung 


Auf Grund der Kenntnis der Eigenwerte von K kommt man durch eine elemen- 
tare Betrachtung, die hier nicht wiedergegeben werden soll, zu der Schranke 
p=2,7. Die Bestimmung von # nach Gleichung (37) ergibt ein nur wenig 
schlechteres Resultat. 
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Berichtigung 

Auf S. 88 der Arbeit von J. H. Hatton “On the efficiency of certain quasi-random 
sequences of points in evaluating multi-dimensional integrals”, erschienen in Band 2, 
S. 84—90 (1960), ist folgendes zu berichtigen: 

(i) Eq. (24): omit the factor 2~*. 

(ii) Line 5: for NA read [NA]. 

(iii) Eq. (27): omit the factor 2—*—, 

(iv) Eq. (28), first line: omit the factor 2—-“-, 
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